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Abstract — Folding a sequence S into a multidimensional box 
is a well-known method which is used as a multidimensional 
coding technique. The operation of folding is generalized in 
a way that the sequence S can be folded into various shapes 
and not just a box. The new definition of folding is based on 
a lattice tiling for the given shape S and a direction in the D- 
dimensional integer grid. Necessary and sufficient conditions 
that a lattice tiling for 5 combined with a direction define 
a folding of a sequence into S are derived. The immediate 
and most impressive application is some new lower bounds 
on the number of dots in two-dimensional synchronization 
patterns. This can be also generalized for multidimensional 
synchronization patterns. The technique and its application 
for two-dimensional synchronization patterns, raise some in- 
teresting problems in discrete geometry. We will also discuss 
these problems. It is also shown how folding can be used to 
construct multidimensional error-correcting codes. Finally, by 
using the new definition of folding, multidimensional pseudo- 
random arrays with various shapes are generated. 

Index Terms — distinct difference configuration, folding, lat- 
tice tiling, pseudo-random array, two-burst-correcting cods 



I. Introduction 

Multidimensional coding in general and two-dimensional 
coding in particular are subjects which attract lot of at- 
tention in the last three decades. One of the main rea- 
sons is their modern applications which have developed 
during these years. Such applications for synchronization 
patterns include radar, sonar, physical alignment, and time- 
position synchronization. For error-correcting codes they 
include two-dimensional magnetic and optical recording as 
well as three-dimensional holographic recording. These are 
the storage devices of the future. Applications for pseudo- 
random arrays include scrambling of two-dimensional data, 
two-dimensional digital watermarking, and structured light 
patterns for imaging systems. Each one of these structures 
(multidimensional synchronization patterns, error-correcting 
array codes, and pseudo-random arrays), and its related cod- 
ing problem, is a generalization of an one-dimensional struc- 
ture. But, although the related theory of the one-dimensional 
case is well developed, the theory for the multidimensional 
case is developed rather slowly. This is due that the fact the 
most of the one-dimensional techniques are not generalized 
easily to higher dimensions. Hence, specific techniques have 
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to be developed for multidimensional coding. One approach 
in multidimensional coding is to take an one-dimensional 
code and to transform it into a multidimensional code. One 
technique in this approach is called folding and it is the 
subject of the current paper. This technique was applied 
previously for two-dimensional synchronization patterns, for 
pseudo-random arrays, and lately for multidimensional error- 
correcting codes. We start with a short introduction to these 
three multidimensional coding problems which motivated 
our interest in the generalization of folding. 
Synchronization patterns 

One-dimensional synchronization patterns were first intro- 
duced by Babcock in connection with radio interference [1]. 
Other applications are discussed in details in [2] and some 
more are given in [3], [4]. The two-dimensional applica- 
tions and related structures were first introduced in [5] and 
discussed in many papers, e.g. [6], [7], [8], [9], [10]. The 
two-dimensional problems has also interest from discrete 
geometry point of view and it was discussed for example 
in [11], [12]. Recent new application in keys predistribution 
for wireless sensor networks [13] led to new related two- 
dimensional problems concerning these patterns which are 
discussed in [14], [15]. It has raised the following discrete 
geometry problem: given a regular polygon with area s on 
the square (or hexagonal) grid, what is the maximum number 
of grid points that can be taken, such that any two lines 
connecting these grid points are different either in their 
length or in their slope. Upper bound technique based on 
an idea of Erdos and Turan [11], [16] is given in [14]. 
Some preliminary lower bounds on the number of dots are 
also given in [14], where the use of folding is applied. 
Folding for such patterns was first used by [10]. An one- 
dimensional ruler was presented as a binary sequence and 
written into a two-dimensional array row by row, one binary 
symbol to each entry of the array. This was generalized 
for higher dimensions, say ni x ri2 x 713 array, by first 
partitioning the array into ri\ two-dimensional arrays of size 
122 x 123. The one-dimensional sequence is written into the 
these ri2 x 77,3 arrays one by one in the order defined by the 
three-dimensional array. To each of these 712 x 77,3 arrays 
the sequence is written row by row. Folding into higher 
dimensions is done similarly and can be defined recursively. 
This technique was used in [10] to generate asymptotically 
optimal high dimensional synchronization patterns. 
Error-correcting codes 

There is no need for introduction to one-dimensional 
error-correcting codes. Two-dimensional and multidimen- 
sional error-correcting codes were discussed by many au- 
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thors, e.g. [17], [18], [19], [20], [21], [22], [23], [24], 
[25]. Multidimensional error-correcting codes are of interest 
when the errors are not random errors. For correction of 
up to t random errors in a multidimensional array, we can 
consider the elements in the array as an one-dimensional 
sequence and use a i-error-correcting code to correct these 
errors. Hence, when we talk about multidimensional error- 
correcting codes we refer to the errors as special ones 
such as the rank of the error array [26], [27], or criss- 
cross patterns [27], [28], [29], etc. An important family of 
multidimensional error-correcting codes are the burst-error- 
correcting codes. In these codes, we assume that the errors 
are contained in a cluster whose size is at most b. The one- 
dimensional case was considered for more than forty years. 
Fire [30] was the first to present a general construction. Op- 
timal burst-correcting codes were considered in [31], [32], 
[33]. Generalizations, especially for two-dimensional codes, 
but also for multidimensional codes were considered in 
various research papers, e.g. [18], [19], [21], [22], [23], [25]. 
In general, "simple" folding of one-dimensional codes were 
not considered for multidimensional error-correcting codes. 
Even so in many of these papers, one-dimensional burst- 
correcting codes and error-correcting codes, were transferred 
into two dimensional codes, e.g. [20], [21], [22], [23], [24], 
[25]. Colorings for two-dimensional coding, which transfer 
one-dimensional codes into multidimensional arrays were 
considered for interleaving schemes [21] and other tech- 
niques [25]. These colorings can be compared to the coloring 
which will be used in the sequel for folding. There is another 
related problem of generating an array in which burst-errors 
can be corrected on an unfolded sequence generated from 
the array [34], [35], [36], [37], [38]. 
Pseudo-random arrays 

The one-dimensional pseudo-random sequences are the 
maximal length linear shift register sequences known as M- 
sequences and also pseudo-noise (PN) sequences [39]. These 
are sequences of length 2 n — 1 generated by a linear feedback 
shift-register of order n. They have many desired properties 
such as 

• Recurrences Property - the entries satisfy a recurrence 
relation of order n. 

• Balanced Property - 2" _1 entries in the sequence are 
ones and — 1 entries in the sequence are zeroes. 

• shift-and-Add Property - when a sequence is added 
bitwise to its cyclic shift another cyclic shift of the 
sequence is obtained. 

• Autocorrelation Property - the out-of -phase value of the 
autocorrelation function is always -1. 

• Window Property - each nonzero n-tuple appears ex- 
actly once in one period of the sequence. 

There are other properties which we will not mention [40]. 
For a comprehensive work on these sequences the reader 
is referred to [39]. Related sequences are the de Bruijn 
sequences of length 2" which are generated by nonlinear 
feedback shift-register of order n. These sequences have the 
window property, i.e., each n-tuple appears exactly once in 
one period of the sequence. 



The two-dimensional generalizations of pseudo-noise and 
de Bruijn sequences are the pseudo-random arrays and 
perfect maps [40], [41], [42], [43], [44], [45]. Pseudo- 
random arrays were also called linear recurring arrays 
having maximum=area matrices by Nomura, Miyakawa, 
Imai, and Fukuda [41] who were the first to construct them. 
Perfect maps and pseudo-random arrays have been used in 
two-dimensional range-finding, in data scrambling, and in 
various kinds of mask configurations. More recently, pseudo- 
random arrays have found other applications in new and 
emerging technological areas. One such application is robust, 
undetectable, digital watermarking of two-dimensional test 
images [46], [47]. Another interesting example is the use 
of pseudo-random arrays in creating structured light, which 
is a new reliable technique for recovering the surface of an 
object. The structured-light technique is based on projecting 
a light pattern and observing the illuminated scene from one 
or more points of view [48], [49], [50], [51]. As mentioned 
in these papers, this technique can be generalized to three di- 
mensions; hence, constructions of three-dimensional perfect 
maps and pseudo-random arrays are also of interest. 

The main goal of this paper is to generalize the 
well-known technique, folding, for generating multidimen- 
sional codes of these types, synchronization patterns, burst- 
correcting codes, and pseudo-random arrays. The general- 
ization will enable to obtain the following results: 

1) Form new two-dimensional codes for these applica- 
tions. 

2) Generalize all the multidimensional codes for any 
number of dimensions in a simple way. 

3) Form some optimal codes not known before. 

4) Make these codes feasible not just for multidimen- 
sional boxes, but also for many other different shapes. 

5) Solve the synchronization pattern problem as a dis- 
crete geometry problem for various two-dimensional 
shapes, and in particular regular polygons. 

It is important to note that folding which was used in other 
places in the literature aim only at one goal. Our folding 
aim is at several goals. Even so, our description of folding 
is simple and very intuitive for all these goals. 

The rest of this paper is organized as follows. In SectionlHl 
we define the basic concepts of folding and lattice tiling. 
Tiling and lattices are basic combinatorial and algebraic 
structures. We will consider only integer lattice tiling. We 
will summarize the important properties of lattices and 
lattice tiling. In Section|III]we will present the generalization 
of folding into multidimensional shapes. All previous known 
folding definitions are special cases of the new definition. 
The new definition involves a lattice tiling and a direction. 
We will prove necessary and sufficient conditions that a 
lattice with a direction define a folding. We first present 
a proof for the two-dimensional case since it is the most 
applicable case. We continue to show the generalization for 
the multidimensional case. For the two-dimensional case the 
proofs are slightly simpler than the slightly different proofs 
for the multidimensional case, we will first consider folding 
in which two consecutive elements in the folded sequence 
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are also adjacent, at least cyclically, in the array. This will 
be generalized to folding in which each two consecutive 
elements in the folded sequence are not necessarily adjacent 
in the array. In Section [IV] we give a short summary on 
synchronization patterns and present basic theorems con- 
cerning the bounds on the number of elements in such 
patterns. In Section [V] we apply the results of the previous 
sections to obtain new type of synchronization patterns 
which are asymptotically either optimal or almost optimal. 
In Section [VT] we discuss folding in the hexagonal grid and 
present construction for synchronization patterns in this grid 
with shapes of hexagons or circles. In Section \VU\ we show 
how folding can be applied to construct multidimensional 
error-correcting codes. In section IVIIII we generalize the 
constructions in [41], [40] to form pseudo-random arrays on 
different multidimensional shapes. Conclusion and problems 
for further research are given in Section |IX] 

II. Folding and Lattice Tiling 

A. Folding 

Folding a rope, a ruler, or any other feasible object 
is a common action in every day life. Folding an one- 
dimensional sequence into a D-dimensional array is very 
similar, but there are a few variants. First, we will summarize 
three variants for folding of an one-dimensional sequence 
sosi • • • s m -i into a two-dimensional array A. The general- 
ization for a D-dimensional array is straightforward while 
the description becomes more clumsy. 

Fl. A is considered as a cyclic array horizontally and 
vertically in such a way that a walk diagonally visits 
all the entries of the array. The elements of the 
sequence are written along the diagonal of the rxi 
array A. This folding works (i.e., all elements of the 
sequence are written into the array) if and only if r 
and t are relatively primes. 

F2. The elements of the sequence are written row by row 
(or column by column) in A. 

F3. The elements of the sequence are written diagonal by 
diagonal in A. 

Example 1: 

Example for Fl: 

Given the M-sequence 000111101011001 of length 15, 
we fold it into a 3 x 5 array with a 2 x 2 window 
property (the extra row and extra column are given for better 
understanding of the folding). 






6 


12 


3 


9 





5 


11 


2 


8 


14 


5 


10 


1 


7 


13 


4 


10 





6 


12 


3 


9 









1 





1 








1 


1 





1 


1 


1 


1 











1 


1 





1 





1 









46 



10 



25 



13 



32 



28 



24 



20 



12 



39 



42 



35 



31 



27 



23 



19 



15 



11 



38 



34 



30 



26 



22 



18 



14 



53 



49 



45 



37 



29 



25 



21 



52 



48 



44 



47 



28 



43 



50 



46 



49 



Fig. 1. Folding by diagonals 



Example for F2: 



The following sequence (ruler) of length 1 3 with five dots 
is folded into a 3 x 5 array 
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Example for F3: 

The following ^-sequence in Z31 : {0,1,4,10,12,17} 
(can be viewed as a cyclic ruler) is folded into an infinite 
array (we demonstrate part of the array with folding into a 
small rectangle is given in bold). Note, that while the folding 
is done we should consider all the integers modulo 31 (see 
Figure [TJ. 

Fl and F2 were used by MacWilliams and Sloane [40] to 
form pseudo-random arrays. F2 was also used by Robin- 
son [10] to fold a one-dimensional ruler into a two- 
dimensional Golomb rectangle. The generalization to higher 
dimensions is straight forward. F3 was used in [14] to obtain 
some synchronization patterns in IP . 

B. Tiling 

Tiling is one of the most basic concepts in combina- 
torics. We say that a D-dimensional shape S tiles the D- 
dimensional space IP if disjoint copies of S cover IP . 

Remark 1: We assume that our shape S is a discrete 
shape, i.e., it consists of discrete points of IP such that 
there is a path between any two points of S which consists 
only from points of S. The shape S in IP is usually not 
represented as a union of points in iP ', but rather as a union 
of units cubes in MP with 2 D vertices in IP . Let A be the 
set of points in the first representation. The set of unit cube 
by the second representation is 

W(h,i2,-,in) '■ ■ ■ ■ Jd) € A} , 
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where 

W(»i,ia,...,iD) = {(*l)ia, ■ ■ • , *d) +£l £1+^262- 

< 6 < 1, 1 < * < £>} , 



and £i is a vector of length D and weight one with a owe in 
the ith position. We omit the case of shapes in M. D which 
are not of interest to our discussion. 

This cover of IP with disjoint copies of S is called tiling 
of IP with S. For each shape S we distinguish one of the 
points of S to be the center of S. Each copy of S in a tiling 
has the center in the same related point. The set T of centers 
in a tiling defines the tiling, and hence the tiling is denoted 
by the pair (T,S), Given a tiling (T, S) and a grid point 
i 2 , . . . , io) we denote by c(i\, ■ ■ ■ , id) the center of 
the copy of S for which (ii,«2j • • • > io) £ S. We will also 
assume that the origin is a center of some copy of iS. 

Remark 2: It is easy to verify that any point of S can 
serve as the center of S. If (T,S) is a tiling then we can 
choose any point of S to serve as a center without affecting 
the fact that (T, S) is a tiling. 

Lemma 1: If (T,S) is a tiling then for any given 
point (ii, 12, . ■ ■ G Z-° the point (ii, i 2 , ■ ■ ■ , id) — 
c(i%, i 2 , ■ ■ ■ ,Id) belongs to the shape S whose center is in 
the origin. 

Proof: Let S\ be the copy of S whose center is 
in the origin and S2 be the copy of S with the point 
(ii,i2) ••• >*£»)• Let (xi, x 2 , ■ ■ ■ , xd) be the point in S\ 
related to the point (i\, i 2 , ■ ■ ■ , «d) in S 2 - By definition, 
(*i,«2, ■■■,«£>) = c(ii,i 2 , ...,«£>) + (xi,x 2 , ■ • ■ ,x D ) and 
the lemma follows. ■ 
One of the most common types of tiling is a lattice tiling. 
A lattice A is a discrete, additive subgroup of the real D- 
space R D . W.l.o.g., we can assume that 



A = {u\Vx + u 2 V2 - 



-u D v D ■ ui,. 



,u D E Z} , (1) 



where {vi, v 2 , ■ • • ,Vd} is a set of linearly independent 
vectors in R D . A lattice A defined by dJ} is a sublattice 
of iP if and only if \y\,x>2, . . . , vd} C 7lP . We will be 
interested solely in sublattices of Z D since our shapes are 
defined in IP . The vectors v\, V2, ■ ■ ■ , Vd are called a base 
for ACZ fl , and the D x D matrix 



G 



W21 



V12 
V22 



vdi vdi 



VlD 
V2D 

VDD 



having these vectors as its rows is said to be a generator 
matrix for A. 

The volume of a lattice A, denoted V(A), is inversely 
proportional to the number of lattice points per unit volume. 
More precisely, V(A) may be defined as the volume of the 
fundamental parallelogram 11(A) in MP, which is given by 



IT(A) d = 



<dv d : < & < 1, , 1 < i < D}. 



There is a simple expression for the volume of A, namely, 
V(A) = |detG|. 



We say that A is a lattice tiling for S if the lattice points 
can be taken as the set T to form a tiling (T,S). In this 
case we have that \S\ = V(A) = \ det G|. 

There is a large variety of literature about tiling and lat- 
tices. We will refer the reader to two of the most interesting 
and comprehensive books [52], [53]. 

Remark 3: Note, that different generator matrices for the 
same lattice will result in different fundamental parallelo- 
grams. This is related to the fact that the same lattice can 
induce a tiling for different shapes with the same volume. A 
fundamental parallelogram is always a shape in M. D which 
is tiled by A (usually this is not a shape in lP and as a 
consequence, most and usually all, of the shapes in lP are 
not fundamental parallelograms). 

Lattices are very fundamental structures in various coding 
problems, e.g. [54], [55], [56] is a small sample which does 
not mean to be representative. They are also applied in 
multidimensional coding, e.g. [21]. This paper exhibits a 
new application of lattices for multidimensional coding and 
for discrete geometry problems. 

Lemma 2: Let A be a D-dimensional lattice, with a 
generator matrix G, and S be a Z?-dimensional shape with 
a point at the origin. A is a lattice tiling for S if and only if 
I det G| = |«S| and there are no two points (ii, 12, ■ . ■ 
and (ji, j 2 , ...,j D )inS such that {ii-Jufa-h, ...,id- 
jd) is a lattice point. 

Proof: Assume first that A is a lattice tiling for S. The 
condition on the volume of S is trivial. Assume the contrary 
that (it, i 3 , ...,id) and (jx,j 2 , ■ ■ ■ ,3d) are in the copy of S, 
whose center is in the origin, and (ii — j\, 12 — j 2 , . . . , id — 
jd) is a lattice point. It follows that the point (i\, i 2 , ■ ■ ■ 
is contained in the shape centered in the origin and also 
in the shape centered at (ii — ji,i 2 — J2> • • • >*d — 3d), a 
contradiction to the fact that A is a lattice tiling for S. 

Now, assume that | det G| = |«S| and there are no two 
points (ii,« 2 , •••,«£>) and (3-1,3-2, ■ ■ ■ ,3d) in S such that 
*2 — 32, ■ ■ ■ , id— 3d) is a lattice point. We choose the 
point of S which is in the origin to be the center of S and we 
place copies of S on each lattice point such that the center 
coincide with the lattice point. Since | det G| = |<S| we only 
have to show that there is no point which is contained in 
two different copies of S in order to complete the proof that 
A is a lattice tiling for S. Assume the contrary that the point 
P is contained in two copies of S with centers at C\ and 
C2. Similarly to the proof of Lemma Q] it can be shown that 
P — C\ and P — C2 are points in the copy of S centered 
at the origin, But, P — Ci - (P — C 2 ) = C 2 — Ci is a 
lattice point (since it is a difference of two lattice points), a 
contradiction to the assumption. Hence, A is a lattice tiling 
for S. M 

Corollary 1: Let A be a £>-dimensional lattice, with a 
generator matrix G, and S be a D-dimensional shape. A is 
a lattice tiling for S if and only if | det G| = |«S| and there 
are no two points (ij, i 2 , ■ ■ ■ , id) and {31,32, ■ ■ ■ ,3d) in any 
copy of S such that (i\ — jx, 12 — 32, ■ ■ ■ ,^d— 3d) is a lattice 
point. 
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III. The Generalized Folding Method 

In this section we will generalize the definition of folding. 
All the previous three definitions (Fl, F2, and F3) are special 
cases of the new definition. The new definition involves 
a lattice tiling A, for a shape S on which the folding is 
performed. 

A ternary vector of length D, (d\, d®, . . . , do), is a word 
of length D, where di S {— 1, 0, +1}. 

Let S be a £>-dimensional shape and let 6 = 
(di, d2, ■ ■ ■ , dn) be a nonzero ternary vector of length D. 
Let A be a lattice tiling for a shape S, and let Si be the 
copy of S which includes the origin. We define recursively a 
folded-row starting in the origin. If the point (ij.,^, ■ • • ,io) 
is the current point of Si in the folded-row, then the next 
point on its folded-row is defined as follows: 

• If the point (ii + d±,i2 + d2, ■ ■ ■ ,iD + dr>) is in Si then 
it is the next point on the folded-row. 

. If the point (ii + d\, i<i + cfe, • • • ,io + do) is in S 2 ^ 
>Si whose center is in the point (ci, C2, . . . , Cd) then 
(ii +di — ci,i 2 + d% — C2, . . . ,in + do — cd) is the 
next point on the folded-row (this is a point in Si by 
Lemma [T). 

The new definition of folding is based on a lattice A, a 
shape S, and a direction S. The triple (A, S, S) defines a 
folding if the definition yields a folded-row which includes 
all the elements of S. It will be proved that only A and 
S determine whether the triple (A, S, S) defines a folding. 
The role of S is only in the order of the elements in the 
folded-row; and of course A must define a lattice tiling for 
S. Different lattice tilings for the same shape S can function 
completely different. Also, not all directions for the same 
lattice tiling of the shape S should define (or not define) a 
folding. 

Remark 4: It is not difficult to see that the three folding 
defined earlier (Fl, F2, and F3) are special cases of the new 
definition. The definition of the generator matrices for the 
three corresponding lattices are left as an exercise to the 
interested reader. 

Remark 5: The definition of ternary vectors for the direc- 
tion, in which the folding is performed, is given to guarantee 
that two consecutive elements in the folded-row, are also 
adjacent (possibly cyclically) in the shape S. 

Example 2: Let S be a 2 x 2 square. Let Ai be the lattice 
whose generator matrix given by the matrix 



Gi = 



2 2 
2 



Ai defines a lattice tiling for S. None of the four possible 
ternary vectors of length 2 define a folding with A (and S). 

Let A2 be the lattice whose generator matrix given by the 
matrix 

" 2 1 

2 



G 2 



A2 also defines a lattice tiling for S. Each one of the 
directions (+1, 0), (+1, +1), and (+1, —1) defines a folding 
with A (and S). Only the direction (0, +1) does not define 
a folding with A (and S). 



How many different folded-rows do we have? In other 
words, how many different folding operations are defined in 
this way? There are 3^ — 1 non-zero ternary vectors. If A 
with the ternary vector (di, d%, . . . , do) define a folding then 
also A with the ternary vector (— di, — c/2, • • • , — do) define a 
folding. The two folded-rows are in reverse order, and hence 
they will be considered to be equivalent. If two folded-rows 
are not equal and not a reverse pair then they will considered 
to be nonequivalent. The question whether for each D, there 
exists a £>-dimensional shape S with 9 different folded- 
rows will be partially answered in the sequel. Meanwhile, 
we present an example for D = 2. 

Before the example we shall define how we fold a 
sequence into a shape S. Let A be a lattice tiling for the 
shape S for which n = \S\. Let 6 be a direction for which 
(A,S,S) defines a folding. Let B = 6061... t> e a 
sequence of length n. The folding of B induced by (A, S, 6) 
is denoted by (A, S, S, B) and defined as the shape S with 
the elements of B, where bi is in the ith entry of the folded- 
row in S defined by (A, S, 6). 

Example 3: Let A be the lattice whose generator matrix 
given by the matrix 



G 



3 2 
7 1 



One can verify that shapes tiled by this lattice have different 
folded-rows. It can be proved that this is the lattice with the 
smallest volume which has this property, i.e., that the four 
folded-rows are different. 

If our shape S is an 1 x 11 array then the folding of a 
sequence with length 1 1 is defined as follows (the position 
labelled with an i is the place of the ith element of the 
sequence). 

For the direction vector (+1,0) the order is given by 
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For the direction vector (+1,-1) the order is given by 
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If our shape S is given by 



then the folding of a sequence of length 1 1 is depicted in 
Figure |2] 

Finally, if our shape S is given by 
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with the direction vector (+1,0); 



with the direction vector (0, +1) 



with the direction vector (+1. +1) 



with the direction vector (+1,-1). 



Fig. 2. Folding of the first shape 
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with the direction vector (+1,0); 



with the direction vector (0, +1); 
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with the direction vector (+1,-1). 



Fig. 3. Folding of the second shape 



then the folding of a sequence of length 1 1 is depicted in 
Figure [3] 

Next, we aim to find sufficient and necessary conditions 
that a triple (A, S, 6) defines a folding. We start with a 
simple characterization for the order of the elements in a 
folded-row. 

Lemma 3: Let A be a lattice tiling for the shape S and 
let S = (dt, d%, . . . , do) be a nonzero ternary vector. Let 
g(i) — (i ■ di, ■ ■ ■ , i ■ do) — c(i ■ dt,..-,i ■ dp) and let 
it, %2 be two integers. Then g(it) — g(ii) if and only if 
g(it + 1) = g(i 2 + 1). 

Proof: The lemma follows immediately from the ob- 
servation that g(ii) = g(ii) if and only if (ii -dt, ■ ■ ■ , it -dr>) 
and (£2 • dt, . . . , i% ■ do) are the same related position in S, 



i.e., corresponds to the same position of the folded-row. ■ 
The next two lemmas are an immediate consequence of 

the definitions and provide a concise condition whether the 

triple (A, S, S) defines a folding. 

Lemma 4: Let A be a lattice tiling for the shape S and let 

5 = (dt,di, . . . , do) be a nonzero ternary vector. (A, S, 6) 
defines a folding if and only if the set {(i • dt, i • di, ■ ■ ■ , i • 
do) — c(i ■ dt,i ■ d2, . . . ,i ■ do) ■ < i < \S\} contains 
|«S| distinct elements. 

Proof: The lemma is an immediate consequence of 
Lemmas [T] [3] and the definition of folding. ■ 
Lemma 5: Let A be a lattice tiling for the shape S and let 

6 = (dt, di, . . . , do) be a nonzero ternary vector. (A, S, 6) 
defines a folding if and only if (|«S| -dt, ■ ■ ■ , |<S| ■ dp) — c(\S\ ■ 
dt,..., \S\ ■ d D ) = (0, • • . , 0) and for each i, < i < \S\ 
we have (i-dt,.-- , i- dx>) — c(i- dt, ■ ■ ■ ,i-do) ^ (0, . . . ,0). 

Proof: Assume first that (A, S, S) defines a folding. If 
for some < j < \S\ we have (j ■ dt, ■■ ■ ,j • do) — c(j ■ 
dt, ■ ■ ■ , j ■ do) — (0, ... ,0) then g(j) — g(0) and hence by 
Lemma [3] the folded-row will have at most j elements of S. 
Since j < \S\ we will have that (A, S, S) does not define 
a folding. On the other hand, Lemma [3] also implies that if 
(A, 5, 5) defines a folding then g(\S\) = (0, . . . , 0). 

Now assume that (|5|-di, . . . , \S\-djj)—c(\S\-dt, . . . ,\S\- 
do) = (0, ■•■,0) and for each i, < i < \S\ we have 
(i ■ dt, ■ ■ ■ , i ■ do) — c(i • dt, . ■ ■ ,i ■ do) 7^ (0, . . . , 0). Let 
< it < i2 < if 9(h) = 9(^2) then by Lemma [3] 
we have g(ii — it) — g(0) — (0, ... ,0), a contradiction. 
Therefore, the folded-row contains all the elements of S 
and hence by definition (A, S, 5) defines a folding. ■ 

Corollary 2: If (A, 5, 6), S = (dt, di, ■ ■ ■ , do), defines a 
folding then the point (|5| -dt, ■ ■ ■ , \S\ -do) is a lattice point. 

Before considering the general 7J-dimensional case we 
want to give a simple condition to check whether the triple 
(A, S, S) defines a folding in the two-dimensional case. For 
each one of the four possible ternary vector we will give 
a necessary and sufficient condition that the triple (A, S, 6) 
defines a folding. 

Lemma 6: Let G be the generator matrix of a lattice A 
and let s = |detG|. Then the points (0, s), (s,0), (s,s), 
and (s, —s) are lattice points. 

Proof: It is sufficient to prove that the points (0, s), 
(s, 0) are lattice points. Let A be a lattice whose generator 
matrix is given by 



G 



Vtl "12 
V21 V 22 

W.l.o.g. we assume that |detG| > 0, i.e., s = vitvn — 
VtiVit. Since U22(«n, ^12) - "12(^21, U22) = (s,0) and 
vtt(v2t,Vn) — vit(vtt,vti) = (0,s), it follows that (0, s), 
(s,0) are lattice points. ■ 
Theorem 1: Let A be a lattice whose generator matrix is 
given by 



G = 



Vtl Vtl 
Vlt V 2 2 

If A defines a lattice tiling for the shape S then the triple 
(A, S, S) defines a folding 
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• with the ternary vector S = (+1,+1) if and only if 
g.c.d.(u 22 - v 2 i,vn - vi 2 ) = 1; 

• with the ternary vector S = (+1,-1) if and only if 
g.c.d.(u 22 + v 2 i,vn + V12) = 1; 

• with the ternary vector S = (+1,0) if and only if 
g.c.d.(wi 2 , V22) = l; 

• with the ternary vector 5 = (0, +1) if and only if 
g.c.d.On, v 2 i) = 1. 

Proof: We will prove the case where S = (+1,+1); 
the other three cases are proved similarly. 

Let A be a lattice tiling for the shape S. By Lemma [6] 
we have that (|«S|, |«S|) is a lattice point. Therefore, there 
exist two integers otx and a 2 such that ai(wn,fi 2 ) + 
a2(i>2i,U22) = (|<S|,|<S|), i.e., aiVn + a 2 v 2 i = a\V\2 + 
012V22 = \S\ = VWV22 — v\2V2\- These equations have 
exactly one solution, a.\ = u 22 — V21 and a 2 = vn — V12. 
By Lemma [5] (A, S, 5) defines a folding if and only if 
= c(|5|,|5|) and for each i, < i < \S\ we 
have (i, i) ^ c(i, i). 

Assume first that g.c.d.(u 22 — V21, vu —V12) = 1. Assume 
for the contrary, that there exist three integers i, /3\, and 
(3 2 , such that /3i(wn,wi 2 ) + /3 2 (^ 2 i, w 22 ) = < 

i < \S\. Hence, flivix + f3 2 v 2 i = $\V\i + P2V22 = i, 
i.e., % = " 11 ^" 12 = 2a. Since a x = v 22 - v 2 i and 
g.c.d.(w 22 — V2i,vn — V12) = 1, it follows that (3i = jai 
and /3 2 = 7a 2 , for some integer 7 > (w.l.o.g. we can 
assume tha 7 > 0). Therefore, i = (3\V\\ + /3 2 u 2 i = 
7ai«ii + 7a 2 w 2 i = 7|5| > |«S|, a contradiction. Thus, it 
follows from Lemma|5]that if g.c.d.(u 22 — w 21 , w 11 — w 12 ) = 1 
then (A, S, S) defines a folding with the ternary vector 
5 =(+!,+!). 

Assume now that (A, S, S) defines a folding with the 
ternary vector S = (+1, +1). Assume for the contrary that 
g.c.d.(w 22 - v 2 x,Vii - U12) = v > 1. Since g.c.d.(w 22 - 
v 2 i,vn - vu) = v > 1, it follows that f3 x = " 22 ~'" 21 
and P2 — Vll ~ Vl2 are integers. Therefore, /3i(un, ^12) + 
/3 2 (v 2 i, V22) — ('TT'i '77') and as a consequence ^ is an 
integer. Hence, by Lemma |5] we have that (A, S, S) does 
not define a folding, a contradiction. Thus, if (A, S, S) 
defines a folding with the ternary vector S = ( + 1, +1) then 

g.C.d.(u22 - U21,Ull - V12) = 1. ■ 

Theorem [T] is generalized for the D-dimensional. This 
generalization will be presented in Theorem [18] given in 
Appendix A. 

There are cases when we can determine immediately 
without going into all the computation, whether (A, S, S) 
defines a folding. It will be a consequence of the following 
lemmas. 

Lemma 7: 

• The number of elements in a folded-row does not 
depend on the point of S chosen to be the center of 
S. 

• The number of elements in a folded-row is a divisor of 
|«S|, i.e., a divisor of V(A). 

Proof: By Lemmas [3] and [5] and the definition of the 
folded-row, if we start the folded-row in the origin then the 
number of elements in the folded-row is the smallest t such 



that t ■ S is a lattice point (since the folded-row starts at 
a lattice point and ends one step before it reaches again a 
lattice point). This implies that the number of elements in a 
folded-row does not depend on the point of S chosen to be 
the center of S. We can make any point of S to be the center 
of S and hence any point can be at the origin. Therefore, all 
folded-rows with the direction 5 have t elements. Any two 
folded-rows are either equal or disjoint. Hence t must be a 
divisor |<S| and t does not depend on which point of S is 
the center. ■ 

The next lemma is an immediate consequence from the 
definition of a folded-row. 

Lemma 8: The number of elements in a folded-row is one 
if and only if S is a lattice point. 

Corollary 3: If the volume of a lattice is a prime number 
then it defines a folding with a direction 6 unless 5 is a 
lattice point. 

By Theorem [18] it is clear that we can determine whether 
the triple (A, S, 5) defines a folding only by the lattice A 
and the ternary direction vector S. The role of S is only 
in the fact that A should be a lattice tiling for S. But, it 
would be easier to examine simpler shapes (like rectangle) 
than more complicated shapes even so they have the same 
lattice tiling A. This leads to an important tool that we will 
use to find an appropriate folding for a shape S' . We will 
use a folding of a simpler shape S with the same volume 
and apply iteratively the following theorem. The proof of the 
theorem is an immediate consequence from the definitions 
of lattice tiling and folding. 

Theorem 2: Let A be a lattice tiling for the D- 
dimensional shape S, let 5 — [d\, d,2, ■ ■ ■ , (Id) be a nonzero 
ternary vector, and (A, S, S) defines a folding. Assume the 
origin is a point in the copy S' of S, i 2 , ■ ■ ■ , id) G S' , 
{i\ + di, 12 + d>2, ■ ■ ■ , %d + do) G S, S' ^ S, and the center 
of S is the point (ci,c 2 , . . . ,cd)- Then A is also a lattice 
tiling for the shape Q = S' U + d\, 12 + d2, . ■ ■ , Id + 
do)} +di - ci,i 2 + d 2 -C2,...,iD+dD- c D )} and 
the triple (A, Q, 8) also defines a folding. 

A. Further generalization of folding 

So far we have used a ternary vector to indicate the 
direction in which the supposed folding is performed. The 
use of a ternary vector is implied by a natural requirement 
that consecutive elements on the folded-row will be also 
consecutive elements in the shape (up to cyclic shift). But, as 
we will see in the sequel, and specifically in the application 
of Sections [IV] and IVIII1 we don't need this requirement. 
This leads for further generalization and modification of 
folding which will yield a better understanding of the 
operation and its properties. 

A direction vector (direction in short) of length D, 
(di, d 2 , • • • , do), is a nonzero word of length D, where 
d{ E Z. The definitions of a folded-row and folding remain 
as before with the exception that instead of a nonzero ternary 
vector we use any nonzero integer direction vector. Also, all 
the results obtained in this section remain true with the same 
proofs. The only exception is Theorem[T]for which we need 
a generalized version which will be given in the sequel. 
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Lemma 9: Let A be a lattice tiling for the shape S. 
Let (di, d2, ■ ■ ■ , do) be a direction vector, ■ ■ ■ >^d) 

be a lattice point, and the point (di, d2, . ■ ■ , do) is in the 
shape S whose center is in the origin. Then the folded-rows 
defined by the directions (di, d%, . . . , do) and (i\ + di, 12 + 
d2, ■ ■ ■ , id + do) are equivalent. 

Proof: Follows immediately from the observation that 
c(ii + di,i 2 +d2,...,id + do) = («i,«2, ■ ■ ■ ^d)- ■ 

In view of Lemma [9] we should examine only the 
\S\ — 1 directions related to the points of S whose cen- 
ter is in the origin. Hence, in the sequel each direction 
5 = (di, g?2, . . . , do) will have the property that the point 
(dx, d2, ■ . ■ , do) will be contained in the copy of S whose 
center is in the origin. One might puzzle how this relates to 
the observation that the necessary and sufficient conditions 
that a direction defines a folding depend only on the gen- 
erator matrix of A and not on 5? The answer is that the 
folded-row itself is defined on the elements of S. Therefore, 
A will have different directions and folded-rows depending 
on the shape S. 

Remark 6: If we consider only the |6>| — 1 directions 
related to the points of S whose center is in the origin, 
some on the ternary direction vectors might not be con- 
sidered (directions which form an equivalent folding will 
be considered). This is another reason for the distinction 
between the definitions of direction vectors (ternary vector 
and integer vector). Each definition has a different purpose. 

Lemma 10: Let A be a lattice tiling for the shape S, 
n = \S\. Let S = (di, c/2, • • • , do) be a direction vector and 
let /0/1 . . . /n-i be its folded-row, where /o = (0, 0, . . . , 0) 
and fi = (di, d2, . . . , dfl). Then the direction 5' = fi 
defines a folding if and only if g.c.d.(i,n) = 1. If the 
direction 6' = fi defines a folding then its folded-row is 
/o/i/2; . . . f n -i, where indices are taken modulo n. 

Proof: By definition and by Lemma [3] we have that 
5' = fi = {i-di,i-d 2 , ■ ■ ■ ,i-dr>) — c(i-d\,i-d2, ■ ■ ■ ,i-d D ) 
and fi.i = (£■ i ■ di,i ■ i ■ d2, ■■■,£■ i ■ do) — c(£ ■ i ■ d\,£- %■ 
d2, ...,£■ i ■ do)- Since the sequence /0/1 . . . f n -i consists 
of n distinct points of 1 D ', it follows that the sequence 
/o/i/2i ■ ■ ■ fn-i consists of 7i distinct points of IP if and 
only if g.c.d.(i, n) = 1. Thus, the lemma follows. ■ 

Corollary 4: Let A be a lattice tiling for the shape S. 
There exists one folding with respect to A if and only if the 
number of nonequivalent folding operations with respect to 
A is — 1^-, where <j)(-) is the Euler function. 

Corollary [4] implies that once we have one folding opera- 
tion with its folded-row, then we can easily find and compute 
all the other folding operations with their folded-rows. It 
also implies that once the necessary and sufficient conditions 
for the existence of one folding in the related theorems 
are satisfied, then the necessary and sufficient conditions 
for the existence of many other folding are also satisfied. 
Nevertheless, Corollary |4] does not guarantee that there will 
be a direction which defines a folding. This fact is shown 
in the next example given in terms of a lemma. 

Lemma 11: Let 7 a positive integer greater than one, a\, 

a2,...,o,£), be nonzero integers, and bi, 62 &d be nonzero 

integers such that either hi — a.- L or bi = aij, for each 



1 < i < D, and \{i : h = ^7, 1 < i < D}\ > 2. Let 
S be a D-dimensional shape and A be a lattice tiling for S 
whose generator matrix is given by 

"61 ... 
b 2 ... 

... b D _ 

Then there is no direction S for which the triple (A, S, 6) 

defines a folding. 

Proof: Let 5 = (d\, di, ■ ■ ■ , dp) be any direction vector 

and let a = 7n£=i a »- Then, a < \S\ and for any given 

shape S for which A is a lattice tiling we have (a ■ d\ , a ■ 

d 2: . . . ,<t- d D ) - c(a ■ di,a ■ d 2 , . . . ,a ■ d D ) = (0, 0, . . . , 0). 

Hence, by Lemma [5] the triple (A, S, 6) does not define a 

folding. ■ 

Lemma 12: Let A be a lattice tiling for the shape S. If \S\ 

\s\—i 

is a prime number then there exists — ^ — different directions 

151 — 1 

which form — L — nonequivalent folded-rows. 

Proof: Let p = \S\ be a prime number. By Corollary [3] 
a direction 6 defines a folding if a and only if S is not a 
lattice point. A shape S in the tiling contains exactly one 
lattice point. Therefore, by Corollary O any one of the p— 1 
directions defined by the non-lattice points of S defines a 
folding. ■ 

Example 4: Consider the lattice A of Example [3] It is a 
lattice tiling for three shapes given in Example [3] For each 
shape, four nonequivalent folding operations are given in 
Example [3] We will demonstrate the fifth one now. 

For the 1 x 11 array the fifth folding operation has the 
direction vector (+2. 0) and the order is given by 
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For the second shape and the direction vector (+2, 0), the 
order is given by 
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For the third shape and the direction vector (+1, +2), the 
order is given by 
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We continue now with the theorem which generalizes 
Theorem[TJ Indeed, it was enough to prove the generalization 
only for the D-dimensional case. But, we feel that making 
the generalizations one step at a time, first for D = 2 and 
after that for any D > 2, will make it easier on the reader, 
and especially as we are using some different reasoning in 
these two generalizations. 

Theorem 3: Let A be a lattice whose generator matrix is 
given by 

q _ vn V12 
V21 v 2 2 
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Let d\ and d 2 be two positive integers and r = 
g.c.d.(<ii, ft). If A defines a lattice tiling for the shape 5 
then the triple (A, 5, 6) defines a folding 

• with the ternary vector 5 = (+di,+d 2 ) if and 
only if g . c A.( dlV22 - d2V21 , d2Vll ~ dlV12 ) = 1 and 
g.c.d.(r, |5|) = 1; 

• with the ternary vector S = (+ft,— ft) if and 
only if g.c.d.( d ' V22 t d2V2 ' , d2Vll + dlVl2 ) = 1 and 
g.c.d.(r, |5|) = 1; 

• with the ternary vector 5 = (+d%,0) if and only if 
g.c.d.(ui2, V22) = 1 and g.c.d.(di, |5|) = 1; 

• with the ternary vector 5 = (0, +d 2 ) if and only if 
g.c.d.(uii, V21) = 1 and g.c.d.(d 2 , |5|) = 1. 

Proof: We will prove the case where 6 — (+di, +^2); 
the other three cases are proved similarly. 

Let A be a lattice tiling for the shape 5. By Lemma [6] 
we have that (|5| • d\, \S\ ■ d 2 ) is a lattice point. Therefore, 
there exist two integers a\ and a 2 such that a\{vu, V12) + 
0t2{v2uV22) = (\S\-di, \S\'d2), i.e., aivu+a 2 V2i = di\S\, 
a\v\2 + OL2V22 — d<x\S\, and |5| = v\\V22 — vi2«2i- These 
equations have exactly one solution, ot\ = d\ 1722 — ^2^21 and 
a 2 = ftvn — d\V\2. By Lemma[5] (A,5, <J) defines a folding 
if and only if (|<S| • di, |5| ■ d 2 ) — c(|5| ■ d\, |5| ■ d 2 ) and for 
each i, < i < |5| we have (t ■ d\, i ■ d 2 ) ^ c(i ■ d\, i ■ d 2 ). 

Assume first that g.c.d.( dlV22 - d2V21 , d ^~ d ^ ) = 

I and g.c.d.(r, |<S|) = 1. Assume for the contrary, 
that there exist three integers i, (3%, and (3%, such that 
/3i(wn,wi 2 ) + ^2(«2i,f22) = (i ■ di,i ■ d 2 ), < i < 
151. Hence we have, ^ = hm=^i2 = 2a. Since 

II Pi a 1 V22 — d,2V2i at 

g.c.d.{ d ' V22 ~ d2V2 ' , d ^~ d ^ ) = 1 it follows that ft = 

1 d 1 v 22 -d 2 v 2 T 1 and ^ = J f d 2 v 11 -d 1 v 12 ^ for some < 7 < T . 

Therefore, we have i ■ di — ftwn + P2V21 = ~" d ]j s \ , i.e., 
i = -^r-- But, since g.c.d.(r, |5|) = 1 it follows that 7 = pr, 
for some integer p > 0, a contradiction to the fact that 
< 7 < t. Hence, our assumption on the existence of 
three integers i, ft, and ft is false. Thus, by Lemma [5] 
we have that if g.c.d. C^ 22 ;'^ 21 , d ^~ d ^ ) = 1 and 
g.c.d.(r, |5|) = 1 then (A, 5,5) defines a folding with the 
direction vector J = {+d\, +^2)- 

Assume now that (A, 5, 5) defines a folding with the 
direction vector 5 = (+di,+d 2 ). Assume for the con- 
trary that g.c.d.( dlV22 ~ d2V2 \ d ^- dlV " ) = Vx > 1 or 
g.c.d.(r, |5|) = V2 > 1. We distinguish now between two 
cases. 

case 1: If g.c.d.( rflt,22 ; d2 " 21 , ^11^12 ) = ^ > 1 then 
ft = d i^-d 2 t> 21 an( j ^ = d2t,11 T -^ 12 are integers. There- 
fore, ft (^1^12) +ft(«2i,«22) T - (^r,^r)- Hence, 
is an integer and for the integers ft = ^"22-^21 and 
ft = d '""- d ^" we have ft> U) « 12 ) + ft> 2 i, i; 22 ) = 

(■^■di, ^^2) and as a consequence by Lemma [5] we have 
that (A,S,S) does not define a folding, a contradiction, 
case 2: If g.c.d.(r, |5|) = ^ 2 > 1 then let ft = d ^~ d ^ 
and ft = d ^~ d ^ . Hence, ft On, Ul2 ) +ft(u 2 i, u 22 ) = 

(~-ft) ^ft)- Clearly, ft, ft, and are integers, and as 
a consequence by Lemma|5]we have that (A, 5, 8) does not 
define a folding, a contradiction. 



Therefore, if (A, 5, 5) defines a folding with the ternary 
vector 5 = (+1,+1) then g.c.d.(w 2 2 - «2i, «ii - V12) = 1- 

■ 

The generalization of Theorem[3]for the _D-dimensional case 
is Theorem [18] given in Appendix A. 

The next lemma is an immediate consequence from the 
definitions on equivalent directions and folded-row. 

Lemma 13: If the directions (ft , <i 2 , . . . , dry) and 
(d'i , d 2 , ■ ■ ■ , d' D ) are equivalent then there exists 
a lattice point i 2 , . . . such that either 

{d' 1 ,d' 2 ,...,d' D ) = (h + di,i 2 + d 2 , . ■ ■ ,id + d D ) or 
(di, d' 2 , . . . ,d' D ) = (ii - di,i 2 - fa,-.., id - do)- 

Lemma 14: Let A be a lattice tiling for a shape 5. If |5| 
is a prime number then there exist 2 ternary direction 
vectors which form folding if and only if there does not exist 
a lattice point (ii, i 2 , . . . , ift, where for each i, 1 < j < D, 
we have < 2. 

Proof: By Lemma |7] if |5| is a prime number, then 
the number of elements in a folded-row for a given ternary 
vector (5 is either one or |5|. By Corollary [3] the number of 
elements is one if and only if 6 is a lattice point. 

If there exist two equivalent directions (ft, . . . , ft>) and 
(d[, . . . ,d' D ) then by Lemma Q~3] we have that (ft — 
d'i, . . . , do — d' D ) is a lattice point, where \di — d[\ < 2 
for each i, 1 < i < D (since |ft < 1 and | < 1). 

If there exists a lattice point (ii, . . . , iu) for which < 
2, 1 < j < D, then there exists two ternary vectors 
(ft,..., do) and (ef^, . . . , d^j) f° r which (ii, . ..,«£>) = 
(dj. .dc-di,). " ■ 

The same result is obtained when |5| is not a prime 
number if the necessary conditions of Theorem [18] are 
satisfied for all the related 2 ternary direction vectors. In 
any case, if there exist a lattice point (ii, i 2 , . . . , where 
for each j, 1 < j < D, we have < 2, then there are 
some related ternary direction vectors which form equivalent 
folding. We can also give an answer to this question by 
finding one ternary direction vector which defines a folding 
and using Corollary [4] 

IV. Bounds on Synchronization Patterns 

Our original motivation for the generalization of the 
folding operation came from the design of two-dimensional 
synchronization patterns. Given a grid (square or hexagonal) 
and a shape 5 on the grid, we would like to find what is 
the largest set A of dots on grid points, |A| = m, located 
in 5, such that the following property hold. All the (™) 
lines between dots in A are distinct either in their length or 
in their slope. Such a shape 5 with dots is called a distinct 
difference configuration (DDC). If 5 is an m x m array with 
exactly one dot in each row and each column than 5 is called 
a Costas array [5]. If 5 is a k x m array with exactly one dot 
in each column then 5 is called a sonar sequence [5]. If 5 is 
a kxn DDC array then 5 is called a Golomb rectangle [7]. 
These patterns have various applications as described in [5]. 
A new application of these patterns to the design of key 
predistribution scheme for wireless sensor networks was 
described lately in [13]. In this application the shape 5 might 
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be a Lee sphere, an hexagon, or a circle, and sometimes 
another regular polygon. This application requires in some 
cases to consider these shapes in the hexagonal grid. F3 was 
used for this application in [14] to form a DDC whose shape 
is a rectangle rotated in 45 degrees in the square grid (see 
Figure [T]). Henceforth, we assume that our grid is Z D , i.e., 
the square grid for D = 2. Since the all the results of the 
previous sections hold for Z?-dimensional shapes we will 
continue to state the results in a D-dimensional language, 
even so the applied part for synchronization patterns is two- 
dimensional. 

We will generalize some of the definition given for DDCs 
in two-dimensional arrays [14] for multidimensional arrays. 
The reason is not just the generalization, but we also need 
these definitions in the sequel. Let A be a (generally infinite) 
D-dimensional array of dots in Z D , and let 771, 772, . . . , T)d 
be positive integers. We say that A is a multi periodic (or 
doubly periodic if D — 2) with period (771, 772, . . . , ?7d) if 
A{i\,i2, ■ ■ ■ ,%d) = A{i\ + r)i,i 2 , ■■■Jd) = A(ii,h + 
772, ■ ■ • , id) = ■ ■ ■ = A{ii,i 2l ■ ■ ■ , id + t/d)- We define the 
density of A to be d/ (HjL^j), where d is the number of 
dots in any 771 x 772 x ••■ x 7773 sub-array of A. Note that 
the period (771 , 772, 7773) might not be unique, but that the 
density of A does not depend on the period we choose. We 
say that a multi periodic array A of dots is a multi periodic 
rix x n 2 x ■ • • n d DDC if every n\ x 112 x • • ■ njj sub-array 
of A is a DDC. 

We write i 2l ■ ■ ■ , id) + S for the shifted copy {(ii + 
71,72 + i' 2 , . . . ,i D + i' D ) ■ («i,«2> ■ ■ ■ ,i' D ) e S} of S. We 
say that a multi periodic array A is a multi periodic S-DDC 
if the dots contained in every shift (ii,i 2 , ■ ■ ■ , id) + S of S 
form a DDC. 

The definition of the density is given based on periodicity 
of a D-dimensional box. If /1 is the density, of the multi 
periodic array A, it implies that given a shape S, the average 
number of dots in any shape S shifted all over A is p,\S\. 
This leads to the following theorem given in [14] for the 
two-dimensional case and which has a similar proof for the 
multidimensional case. 

Theorem 4: Let S be a shape, and let A be a multi 
periodic S-DDC of density ix. Then there exists a set of 
at least |~/x|S|] dots contained in S that form a DDC. 

Another important observation from the definition of multi 
periodic S-DDC is the following lemma from [14]. 

Lemma 15: Let A be a multi periodic S-DDC, and let 
S' C S. Then A is a multi periodic S'-DDC. 

Let Si, £>2, ■ ■ ■ be an infinite sequence of similar shapes 
such that |S;+i| > \Si\. Using the technique of Erdos and 
Turan [11], [16], for which a detailed proof is given in [14], 
one can prove that 

Theorem 5: An upper bound on the number of dots in Si, 
7 — »• 00, is lirn,_ ) . 0o (- N /|S7| + o(y/\ST\)). 

Let S and S' be two-dimensional shapes in the grid. We 
will denote by A(S,S') the largest intersection between S 
and S' in any orientation. Our bounds on the number of dots 
in a DDC with a given shape are based on the following 



result. 

Theorem 6: Assume we are given a multi periodic S- 
DDC array A with density /i. Let Q be another shape on 
IP . Then there exists a copy of Q on 1 D with at least 
[> ■ A(S, Q)] dots. 

Proof: Let Q' be the shape such that S'=5nQ and 
\Q'\ = A(S, Q). By Lemma [TBI we have that A is a multi 
periodic Q'-DDC. By Theorem [4] there exists a set of at 
least [~//|Q'|] dots contained in S that form a DDC. Thus, 
there exists a copy of Q on Z D with at least \/i ■ A(S, Q)~\ 
dots. ■ 

In order to apply Theorem [6] we will use folding of 
sequences defined as follows. Let A be an abelian group, and 
let B = {bi, b 2 , . . . , b m } C A be a sequence of m distinct 
elements of A. We say that B is a B 2 -sequence over A if all 
the sums ai 1 + a,i 2 with 1 < i\ < i 2 < m are distinct. For a 
survey on £?2-sequences and their generalizations the reader 
is referred to [57]. The following lemma is well known and 
can be readily verified. 

Lemma 16: A subset B = {a\, a 2 , ... , a m } C A is a B 2 - 
sequence over A if and only if all the differences — a,i 2 
with 1 < 7i ^ i 2 < m are distinct in A. 

Note that if B is a ^-sequence over Z„ and a € Z n , 
then so is the shift a + B = {a + e : e £ B}. The following 
theorem, due to Bose [58], shows that large ^-sequences 
over Z„ exist for many values of n. 

Theorem 7: Let q be a prime power. Then there exists a 
£?2-sequence ai,a 2 , . . . , a m over Z n where n = q 2 — 1 and 

777 = q. 

A. A Lattice Coloring for a Given Shape 

In this subsection we will describe how we apply folding 
to obtain a DDC with a shape S and a multi periodic 
S-DDC. Let A be a lattice tiling for S and let 6 = 
(di, d 2 , . . . , dr>) be a direction vector such that (A,S, S) 
defines a folding. We assign an integer from Z„, n = \S\, to 
each point of Z-° . The lattice coloring C(A, 8) is defined as 
follows. We assign to the point (0,0,..., 0) and we color 
the next element on the folded-row with 1 and so on until 
|S| — 1 to the last element on the folded-row. This complete 
the coloring of the points of the shape S whose center is 
the origin. To position (71,72, . . . ,id) we assign the color 
of position i 2 , . . . , ijj) — c(ii,i 2 , . . . , 773). The color of 
position %2, ■ ■ ■ , Id) will be denoted by C(ii, i 2 , ■ ■ ■ , io)- 

We will generalize the definition of folding a sequence 
into a shape S by the direction S, given the lattice tiling 
A for S. The folding of a sequence B — 60^1 ■ ■ ■ &n-i 
into an array colored by the elements of Z n is defined by 
assigning the value 6j to all the points of the array colored 
with the color 7. If the coloring was defined by the use 
of the folding as described in this subsection, we say that 
the array is defined by (A,S, 8, B). Note, that we use the 
same notation for folding the sequence B into the shape S. 
The one to which we refer should be understood from the 
context. 

Given a point 12, . . ■ , id) € Z , we say that the set of 
points {(ii+£-di, i 2 +£-d 2 , iD+i-do) ■ £ G Z} is a row 
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oflP defined by 8. This is also the row of (it, £ 2 , ■ ■ ■ , £d) 
defined by 8 = (dt, <f 2 , ■ ■ ■ , do)- 

Lemma 17: If the triple (A, <S, 8) defines a folding then 
in any row of TP defined by 8 there are lattice points. 

Proof: Given a point (it, £ 2 , . ■ ■ , Id) and its color 
C(it,i2, ■ ■ ■ , id), then by the definitions of the folding and 
the coloring we have that C(ii + di,i 2 + d 2 , iD + du) = 
C(it,i 2 , ■ ■ ■ , in) + 1 (mod |<S|). Hence, the row defined by 
S has all the values between and |<S| — 1 in their natural 
order modulo |«S|. Therefore, any row defined by 6 has lattice 
points (which are exactly the points of this row which are 
colored with zeroes). ■ 

Corollary 5: If (it, £ 2 , . . . , i D ), (it +et, £ 2 + e 2 , ■ ■ • , id + 
eu), (jt, 32, ■ ■ -,3d), and (jt + e x ,32 + e 2, ■ ■ ■ , 3d + e D ) 
are four points of T D then C(i\ + et, £2 + e 2 , ■ ■ ■ > £d + 
eo)-C(h,i2, • ••,£d) = C(ji+ei, j 2 + e 2 , . . . Jd + cd)- 
C(ji,h, ■ ■ -,3d) (mod \S\). 

Proof: By Lemma [TTI to each one of these four points 
there exists a lattice point in its row defined by 8. Let 

• Pi = (it + a x ■ di,i 2 + at ■ d2, . ■ ■ ,tD + o>t ■ do) be 
the lattice point in the row of (£1,12, • • • , «r>); 

• P2 = (ji + a 2 ■ dt, 32 + a 2 ■ d 2 , . . . , 3d + a 2 ■ d D ) the 
lattice point in the row of (31,32, ■ ■ ■ ,3d)', 

• P. 3 = ((ii+ei) + a 3 -di,(i 2 + e 2 ) + a 3 -d 2 ,...,(i D + 
eo) + a 3 • do) the lattice point in the row of (it + 
et,i 2 + e 2 , • ■ • , id + e D ). 

Therefore, P 4 = P 2 + P 3 - Pt = ((ji + ex) + (a 2 + 
a 3 - at) ■ dt, (32 + e 2 ) + (a 2 + a 3 - at) ■ d 2 , . . . , (jo + 
en) + (a 2 + a 3 — at) ■ c?_d) is also a lattice point. P4 
is a lattice point in the row, defined by 8, of (ji + 
ei, 32 + e 2 , . . . , 3d + e_o)- All these four points are colored 
with zeroes. Hence, C(it, £ 2 , . . . , io) = —at (mod |<S|), 
C(h + ei,£ 2 + e 2 ,...,£z3 + e D ) = -a 3 (mod 
C(Ju32,---Jd) = ~a% (mod and C(jt + ei,j 2 + 
e 2 , . . . ,30 + eo) = -(a 2 +a 3 - a x ) (mod |<S|). Now, the 
claim of the corollary is readily verified. ■ 
Corollary 6: If 8' is an integer vector of length D then 
there exists an integer e(8') such that for any given point 
P = {h,i 2 ,...,i D ) we have C(P + 8') = C(P) + 
e(8') (mod 

Corollary 7: If the triple (A, S, 8) defines a folding and 
B is a _B 2 -sequence over Z n , where n = \S\, then the array 
A defined by (A, S, 8, B) is multi periodic. 

Proof: Clearly, the array has period (|«S|,|«S|,...,|«S|) 
and the result follows. ■ 

Theorem 8: If the triple (A, S, 8) defines a folding and B 
is a £? 2 -sequence over Z„, where n = \S\, then the pattern 
of dots defined by (A, S, 8, B) is a multi periodic 5-DDC. 

Proof: By Corollary [7] the constructed array is multi 
periodic. 

Since (A, <S, 8) defines a folding it follows that the |«S| 
colors inside the shape S centered at the origin are all dis- 
tinct. By Corollary |5J for the four positions (£1, £ 2 , . . . , £d), 
(it + et,i 2 + e 2 ,...,i D + e D ), (31,32, ■ ■ - ,3 d), and (ji + 
ei,j2 + e 2 ,...,j D + e D ) we have that C(i x + et,i 2 + 
e 2 ,.-.,iD + e D ) - C(£i,£ 2 , . . . ,i D ) = C(ji + et,j 2 + 
e 2 , • Jd + sd) - C(jt,j 2 , . . .,j D ) (mod |5|). Hence, at 



most three of these integers (colors) are contained in B. It 
implies that if these four points belong to the same copy of 
S on the grid then at most three of these points have dots, 
since the dots are distributed by the B 2 -sequence B. Thus, 
any shape S on IP will define a DDC and the theorem 
follows. ■ 

Corollary 8: If the triple (A, S, 8) defines a folding and 
B is a £? 2 -sequence over Z„, where n = |<S|, then the pattern 
of dots defined by (A, S, 8, B) is a DDC. 
Note, that the difference between Theorem[8]and Corollary [8] 
is related to the folding into TL D and S, respectively. The last 
lemma is given for completeness. 

Lemma 18: If (A, S,8) defines a folding then the |<S| 
colors inside any copy of S on a lP are all distinct. 

Proof: Let S x and S 2 be two distinct copies of S 
on IP . Clearly, S 2 = (et, ■ ■ ■ ,eo) + St- By Corollary [5] 
for each (it, . . . , i D ), (ji, . . . , j D ), e Si we have C(it + 
et,...,i D + e D ) - C(it, .. . ,i D ) = C(ji + e x , . . . , j D + 
en) — C(jt, . . . , 3d) (mod |«S|). Therefore, if Si contains 
|<S| distinct colors then also S 2 contains |<S| distinct colors. 
The lemma follows now from the fact that (A, S, 8) defines 
a folding and therefore all the colors in the shape S whose 
center is in the origin are distinct. ■ 
Note, that theorem [8] is also an immediate consequence of 
Lemma [18] 

V. Bounds for Specific Shapes 

In this section we will present some lower bounds on the 
number of dots in some two-dimensional DDCs with specific 
shapes. In the sequel we will use Theorem[6] Theorem[S] and 
Corollary [8] to form DDCs with various given shapes with a 
large number of dots. To examine how good are our lower 
bounds on the number of dots, in a DDC whose shape is 
Q, we should know what is the upper bound on the number 
of dots in a DDC whose shape is Q. By Theorem [5] we 
have that for a DDC whose shape is a regular polygon or 
a circle, an upper bound on the number of dots is at most 
V / s + o(v / s), where the shape contains s points of the square 
grid and s — * 00. One of the main keys of our constructions, 
and the usage of the given theory, is the ability to produce 
a multi periodic 5-DDC, where S is a rectangle, the ratio 
between its sides is close to any given number 7, and if 
its area is s then the number of dots in it is approximately 
>/S + o(VS). 

For the construction we will need the well known Dirich- 
let's Theorem [59, p. 27]. 

Theorem 9: If a and b are two positive relatively primes 
integers then the arithmetic progression of terms ai + b, for 
i = 1,2, contains an infinite number of primes. 

The following theorem is a well known consequence of 
the well known Euclidian algorithm [59, p. 11]. 

Theorem 10: If a and f3 are two integers such that 
g.c.d.(a, 0) = 1 then there exist two integers c a and cp 
such that c a a + cp(3 = 1. 

The next theorem makes usage of these well known old 
foundations. 

Theorem 11: For each positive number 7 and any e > 0, 
there exist two integers nt and n 2 such that 7 < ^1 < 7 + e; 
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and there exists a multi periodic 5-DDC with \J a ■ bR+o(R) 
dots whose shape is an n% x n-i = (aR+o(R)) x (bR+o(R)) 
rectangle, where n\n-i — p 2 — 1 for some prime p, and n\ 
is an even integer. 

Proof: Given a positive number 7 and an e > 0, it is 
easy to verify that there exist two integers a and (3 such that 
V7 < I < V7+~e and g.c.d.(a,/3) = 2. By Theorem [10] 
there exist two integers c a , cr such that either c a a + 2 = 
cp(3 > or c/3/3 + 2 = c Q o; > 0. 

Assume c a a + 2 — cp[3 > (the case where + 
2 = c a a > is handled similarly). Clearly, any factor of a 
cannot divide c a a+l. Since j3 divides c Q o+2, it follows that 
a factor of (3 cannot divide c a a+ 1. Hence, g.c.d.(a/3, c a a + 
1) = 1. Therefore, by Theorem [9] there exist infinitely many 
primes in the sequence aj3R + c a a + 1, R = 1, 2, . . . . 

Let p be a prime number of the form a(3R+c a a+l. Now, 
p 2 -l = = (a/3i? + c Q a + 2)(a/3i? + c Q o;) = 

(a/3i? + cp/3){a/3R + c a a) = {a 2 R + acp)(f3 2 R + /3c a ). 
Thus, a (f3 2 R + f3c a ) x (a 2 i? + 00/3) rectangle satisfies the 
size requirements for the n\ x 712 rectangle of the Theorem. 

Let a = /3 2 , b = a 2 , n\ = (3 2 R + (3c a , n 2 = a 2 R + acp, 
and let S be an n\ x n 2 rectangle. Let A be the a lattice 
tiling for S with the generator matrix 



G 



n-2 




"1 
2 



where 6 = 1 if n\ = (mod4) and 6* = 2 if 77,1 = 2 (mod 
4). By Theorem^ (A,S,5), 8 = (+1,0), defines a folding. 

The existence of a multi periodic 6>-DDC with V a • fei? + 
o(i?) dots follows now from Theorems [7] and |8] ■ 

The next key structure in our constructions is a certain 
family of hexagons defined next. A centroid hexagon is an 
hexagon with three disjoint pairs of parallel sides. If the four 
angles of two parallel sides (called the bases of the hexagon) 
are equal and the four other sides are equal, the hexagon 
will be called a quasi-regular hexagon and will be denoted 
by QRH(kj, b, h), where b is the length of a base, h is the 
distance between the two bases, and b + 2w is the length 
between the two vertices not on the bases. We will call the 
line which connects these two vertices, the diameter of the 
hexagon (even if it might not be the longest line between 
two points of the hexagon). Quasi-regular hexagon will be 
the shape S that will have the role of S when we will apply 
Theorem [6] to obtain a lower bound on the number of dots 
in a shape Q which usually will be a regular polygon. In 
the sequel we will say that — « 7, when we means that 

7 < I < 7 + e- 

We want to show that there exists a quasi-regular 
hexagon QRH(w, b, h) with approximately y/ (b + w)h + 
(V (P + w)h) dots. By Theorem [TT] there exists a doubly 
periodic 5-DCC, where S is an n\ x n 2 = (aR + o(R)) x 
(f3R + o{R)) rectangle, such that ^ » ^±H, ni n 2 =p 2 -l 
for some prime p, and n\ is an even integer. The lattice 
A of Theorem [TT] is also a lattice tiling for a a shape S', 
where S' is "almost" a a quasi-regular hexagon QRH(u>, b, h) 
(part of this lattice tiling is depicted in Figure 0). By 
Theorem [TJ (A, S, 5), 5 = (+1,0), defines a folding for 
this shape too. Hence, we obtain a doubly periodic 6>'-DCC, 









/ 


















Fig. 4. From rectangle to "almost" quasi-perfect hexagon with the same 
lattice tiling 



where S' is "almost" a a quasi-regular hexagon QRH(w, b, h) 
with approximately y/{b + w)h + o(\J (b + w)h) dots. This 
construction implies the following theorem. 

Theorem 12: A lower bound on the number of dots in 
a regular hexagon with sides of length R is approximately 



V2 



-R + o{R). 



Now, we can give a few examples for other specific 
shapes, mostly, regular polygons. To have some comparison 
between the bounds for various shapes we will assume that 
the radius of the circle or the regular polygons is R (the 
radius is the distance from the center of the regular polygon 
to any one its vertices). We also define the packing ratio 
as the ratio between the lower and the upper bounds on the 
number of dots. The shape S that we use will always by a 
multi periodic 5-DDC on a multi periodic array A. 

A. Circle 

We apply Theorem [6] where S is a regular hexagon 
with radius p and Q is a circle with radius R, sharing 
the same center. The upper bound on the number of dots 
in Q is ^/ttR + o(R). A lower bound on the number of 

/3V3 



dots in S is approximately 
density of A is approximately f 2 



V2 



o(p) and hence the 
. Let 9 be the angle 



between two radius lines to the two intersection points of 
the hexagon and the circle on one edge of the hexagon. 
We have that A(5, Q) = (vr - 30 + 3 sin 0)R 2 and p = 

cos — 

i-R. Thus, a lower bound on the number of dots in Q 



-A(<S, Q). The maximum is obtained when 



V2|5| 

9 = 0.536267 yielding a lower bound of 1.70813i? + o(R) 
on the number of dots in Q and a packing ratio of 0.9637. 

We must note again, that even so this construction works 
for infinitely many values of R, the density of these values 
is quite low. This is a consequence of Theorem [TT] which 
can be applied for an arbitrary ratio 7 only when the 
corresponding integers obtained by Dirichlet's Theorem are 
primes. Of course, there are many possible ratios between 
the sides of the rectangle that can be obtained for infinitely 
many values. A simple example is for any factorization of 
p 2 — 1 = ?7.i7i2 we can form an n\ x rt 2 DDC and from 
its related quasi-regular hexagons. We won't go into details 
to obtain bounds which hold asymptotically for any given 
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R as we conjecture that the construction for quasi-regular 
hexagon can be strengthen asymptotically for almost all 
parameters. Nevertheless, we will show briefly how we can 
use a doubly periodic 6>-DDC, where S is a square to obtain 
a lower bound for the number of dots in a DDC whose shape 
is a circle. We use a doubly periodic S-DDC, where S is a 
(p + 1) X (p — 1) rectangle. For a lattice tiling of S we use 
a lattice A with the generator matrix 

r p-i H+i + e - 

[ p+l 

where 9 = 1 if p + 1 = (mod 4) and 9 = 2 if p + 1 = 
2 (mod 4). By Theorem Q] (A,S,S), S = (+1,0), defines 
a folding. We can use Theorem [2] to obtain a new shape 
S' which produces better intersection with a circle, and a 
better lower bound on the number of dots in it (the previous 
best packing ratio obtained with the method implied only 
by Theorem [6] (without using Theorem [2] and better multi 
periodic 5-DDCs) was 0.91167 and it was given in [14]). 



B. Regular Polygon 

For regular polygons with small number of sides we will 
use specific constructions which are given in Appendix C. 
If the number of sides is large we will use Theorem [6] 
where Q will be the regular polygon and S is a regular 
hexagon. Assume that the regular polygon has n sides, R is 
its radius, and p is the radius of the regular hexagon. The 
area of the hexagon is ^^-p 2 and hence the density of the 
doubly periodic array A is approximately f 2 . The area 

V 3\/3p 

of the regular polygon is " R " n " and hence an upper 

bound on the number of dots in Q is " R + o{R). 

For simplicity we will further assume that n = 12k (the 
results for other values of n are similar, but the constructions 
become slightly more complicated for short description. We 
will choose a regular hexagon which has a joint center with 
the regular polygon. We further choose it in a way that S and 
Q intersect in exactly 12 vertices of Q equally spread. We 
will also make sure that each side of S intersects exactly 
two vertices of Q with equal distance from the nearest 
vertices of S to these two intersection points. It implies 
that A(6>, Q) = 6+n " n " R 2 and hence a lower bound on 
the number of dots is 6+ " hm a j{ -|_ o(R). Some values 

2-34(v/3+l) 

obtained from this construction are given in Table U 

For small values of n, specific constructions are given in 
Appendix C. For some constructions we need DDCs with 
other shapes like a Corner and a Flipped T which are defined 
in Appendix B, where also constructions of multi periodic 
5-DDCs for these shapes are given. Table J] summarizes the 
bounds we obtained for regular polygons and a circle in the 
square grid. The same techniques can be used for any D- 
dimensional shape. Finally, we note that the problem is of 
interest also from discrete geometry point of view. Some 
similar questions can be found in [12]. 



TABLE I 

Bounds on the number of dots in an n-GON DDC 



n 


upper bound 


lower bound 


packing ratio 


3 


1 


13975.R 


1.02462/? 


0.899 


4 


1 


41421R 


1.41421i? 


1 


5 


1 


54196R 


1.45992R 


0.9468 


6 


1 


61185R 


a 1.61185R 


Si 1 


7 


1 


65421R 


1.58844i? 


0.960241 


8 


1 


68179R 


1.62625R 


0.966977 


9 


1 


70075.R 


1.63672R 


0.96235 


10 


1 


71433R 


1.64786i? 


0.961229 


12 


1 


73205/J 


1.66871i? 


0.963433 


24 


1 


76234i?, 


1.69815R 


0.963578 


36 


1 


76796/J 


1.70367R 


0.963636 


48 


1 


76992_R 


1.7056iJ 


0.963658 


60 


1 


77083.R 


1.7065.R 


0.963669 


72 


1 


77133R 


1.70699i? 


0.963675 


84 


1 


77163R 


1.70728R 


0.963679 


96 


1 


77182R 


1.70747H 


0.963681 


circle 


1 


77245iJ 


1.70813i? 


0.963708 



VI. Folding in the Hexagonal Grid 

The questions concerning DDCs can be asked in the 
hexagonal grid in the same way that they are asked in 
the square grid. Similarly, they can be asked in dense D- 
dimensional lattices. In this section we will consider some 
part of our discussion related to the hexagonal grid. The 
hexagonal grid is a two-dimensional grid and hence we will 
compare it to Z 2 . In Z 2 there are four different ternary 
direction vectors, while in the hexagonal grid there are three 
different related directions. But, the total number of direc- 
tions depend on the shape in both grids (see Subsection llll-AI 
and especially Corollary [4j. We can define a folded-row and 
folding in the hexagonal grid in the same way as they are 
defined in Z 2 . To prove that the results remain unchanged 
we will describe the well known transformation between the 
hexagonal grid and Z 2 . 

The hexagonal grid is defined as follows. We start by 
tiling the plane R 2 with regular hexagons whose sides have 
length 1 / \/3 (so that the centers of hexagons that share an 
edge are at distance 1). The center points of the hexagons 
are the points of the grid. The hexagons tile R 2 in a way 
that each point (i, 0), i S Z, is a center of some hexagon. 

The transformation uses an isomorphic representation of 
the hexagonal grid. Each point (x, y) 6 Z 2 has the following 
neighboring vertices, 

{(x + a,y + b) | a,be {-1, 0, 1}, a + b ^ 0}. 

It may be shown that the two representations are isomorphic 
by using the mapping £ : R 2 — > R 2 , which is defined by 
£(x, y) = (x + -^,- 2 ^). The effect of the mapping on 
the neighbor set is shown in Fig. [5] From now on, slightly 
changing notation, we will also refer to this representation 
as the hexagonal grid. Using this new representation the 
neighbors of point are 

{(i -1,3- 1), (i - 1, j), (», j ~ 1), (», j + 1), 
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Fig. 6. From a comer CR(9, 9; 5, 4) to hexagonal sphere with radius 4 



Fig. 5. The hexagonal model translation 

Lemma 19: Two lines differ in length or slope in one 
representation if and only if they differ in length or slope in 
the other representation. 

Proof: This claim can be verified easily by observing 
that two lines are equal in length and slope in one represen- 
tation if and only if they are equal in length and slope in 
the other representation. ■ 

Corollary 9: A shape S is a DDC in the hexagonal grid 
if and only if = {£(p) : p e S} is a DDC in Z 2 . 

Clearly, the representation of the hexagonal grid in terms 
of Z 2 implies that all the results on folding in the square 
grid hold also in the hexagonal grid. We will consider 
now the most important families of DDCs in the hexagonal 
grid, regular hexagons and circles. A regular hexagon in the 
hexagonal grid is also called an hexagonal sphere with radius 
R. It is a shape with a center hexagon which includes all 
the points in the hexagonal grid which are within Manhattan 
distance R from the center point. Applying the transforma- 
tion £ on this sphere we obtain a new shape in the square 
grid. This shape is a (2R+1) x (2R+ 1) square from which 
isosceles right triangle with sides of length R are removed 
from the left upper corner and the right lower corner. For the 
construction we use as our shape S, in Theorem [6j a corner 
CR(2i?, w\ +w 2 ; R, w^), where ^- rs 1, \w\ — w%\ < 3 and 
g.c.d.(wi, W2) = 1. In Appendix B a construction for doubly 
periodic 5-DDC, where S is such corner, is given where the 
number of dots in S is approximately a/|«S| + o(yf\S\). By 
Theorem |2] the lattice tiling for S is also a lattice tiling 
for the shape S' obtained from S by removing an isosceles 
right triangle with sides of length R from the lower left 
corner and adding it to the upper right corner of the S 
(see Figure [6)l. The constructed doubly periodic S' -DDC 
can be rotated by 90 degrees or flipped either horizontally 
or vertically to obtain a doubly periodic Q-DDC, where Q 
is approximately an hexagonal sphere with radius R. This 
yields a packing ratio approximately 1 between the lower 
bound and the upper bound on the number of dots. Now, 
it is easy to verify that the same construction, for a DDC 
with a circle shape, given in Subsection IV-AI for the square 
grid will work in the hexagonal grid. For this construction 
we will use regular hexagon and a circle in the hexagonal 
grid to obtain a packing ratio between the lower bound and 
the upper bound on the number of dots in the circle which 
is the same as in the square grid. 

VII. Application for Error-Correction 

In this section we will discuss the usage of folding 
to design optimal (or "almost" optimal) codes which can 



correct adjacent errors in a multidimensional array, i.e., 
a multidimensional 2-burst-correcting code. The construc- 
tion is a generalization of the construction of optimal 
one-dimensional 2-burst-correcting codes given by Abram- 
son [31]. His construction was generalized for larger bursts 
by [32] and [33] who gave a comprehensive treatment 
for this topic. Multidimensional generalization for the 2- 
burst-correcting codes were given in [23], [60]. We will 
give a multidimensional generalization only for the 2-burst- 
correcting codes. The parity-check matrix of a code of length 
2 m _ 1 and redundancy rn + 1, consists of the 2 m — 1 
consecutive nonzero elements (powers of a primitive element 
a) of GF(2 m ) followed by a row of ones. The received 
word has one or two errors depending if the last entry of its 
syndrome is one or zero, respectively. The position of the 
error is determined by the first m entries of the syndrome. 

The generalization of this idea is done by folding the 
nonzero elements of GF(2™ 1 ) into the parity-check matrix 
of a multidimensional code row by row, dimension by 
dimension. Assume that we have a Z?-dimensional array of 
size n\ x ni x ■ • • x nr> and we wish to correct any D- 
dimensional burst of length 2 (at most two adjacent positions 
are in error). The following construction given in [60] is 
based on folding the nonzero elements of a Galois field with 
characteristic 2 into a parity check matrix, where the order 
of the elements of the field is determined by a primitive 
element of the field. 

Construction A: Let a be a primitive element in GF(2 m ) 
for the smallest integer m such that 2 m — 1 > Ylg=i n £- 
Let d = |~log 2 D] and i = 12, . . ■ , id), where < ig < 
ni — 1. Let A be a d x D matrix containing distinct binary 
(i-tuples as columns. We construct the following n\ x n 2 x 
• • • X nry X (m + d + 1) parity check matrix H. 



1 

Ai T mod 2 



for all i = (ii,i2, ■ ■ ■ , Id), where < it < ni — 1. 

The following two theorems were given in [60]. 

Theorem 13: The code constructed in Construction A can 
correct any 2-burst in an n\ x n 2 x • • • x no array codeword. 

Theorem 14: The code constructed by Construction A has 
redundancy which is greater by at most one from the trivial 
lower bound on the redundancy. 

The same construction will work if instead of a D- 
dimensional array our codewords will have have a shape 
S of size 2 m — 1, there is a lattice tiling A for S, and there 
is a direction vector 5 such that (A, S, S) defines a folding. 
The nonzero elements of GF(2 m ) will be ordered along the 
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folded-row of S. Since usually the number of elements in S 
is not 2 m — 1 we should find a shape S' which contains S 
and |<S'| = 2 m — 1. We design a code with the shape of S' 
and since S C S' the code will be able to correct the same 
type of errors in S. 

Finally, the construction can be generalized in a way that 
the multidimensional code will be able to correct other types 
of two errors in a multidimensional array [60]. 

VIII. Application for Pseudo-Random Arrays 

MacWilliams and Sloane [40] gave the name pseudo- 
random sequence to a maximal length sequence obtained 
from a linear feedback shift register. These sequences called 
also PN Pseudo Noise sequences or M-sequences have 
many desired properties as described in [39], [40]. The 
term pseudo-random array was given by MacWilliams and 
Sloane [40] to a rectangular array obtained by folding a 
pseudo-random sequence S into its entries. The constructed 
arrays can be obtained also as what is called maximum- 
area matrices [41]. In [40] it was proved that if a pseudo- 
random sequence of length n = 2 klk2 — 1 is folded into an 
n\ x n2 array such that n\ = 2 kl — 1 > 1, n.2 = ^- > 
1, and g.c.d(?ii, 7i 2 ) > 1 then the constructed array has 
many desired properties and hence they called this array 
A a pseudo- random array. Some of the properties they 
mentioned are as follows: 

1) Recurrences - the entries satisfy a recurrence relation 
along the folding. 

2) Balanced - 2 klk2 ~ 1 entries in the array are ones and 
2fcifc 2 -i _ i entries in the array are zeroes. 

3) Shift-and-Add - the sum of A with any of its cyclic 
shifts is another cyclic shift of A. 

4) Autocorrelation Function - has two values: n in-phase 
and -1 out-of -phase. 

5) Window property - each of the 2 klk2 — 1 nonzero 
matrices of size fci x k% is seen exactly once as a 
window in the array. 

All these properties except for the window property are a 
consequence of the fact that the elements in the folded-row 
are consecutive elements of an M-sequence S. Before we 
examine whether an array of any shape, obtained by folding 
S into it, has these properties we have to define what is a 
cyclic shift of any given shape S (even so we used the term 
without definition before). Our definition will assume again 
that there exists a lattice tiling A for S and a direction S 
such that (A, S, 8) defines a folding. A cyclic shift of the 
shape S (placed on the grid) is obtained by taking the set 
of elements {x + 8 : x G S}. 

Lemma 20: The shape of a cyclic shift of S is S. 

Proof: The cyclic shift is just a shift by 8 of S on the 
grid. Therefore, the shape obtained is also S. ■ 

Theorem 15: Let A be a lattice tiling for a shape S and 
let 8 be a direction such that (A, S, 8) defines a folding. 
If an M-sequence S is folded into S in the direction 8 
then the Recurrences, Balanced, Shift-and-Add, and the 
Autocorrelation Function properties hold for the constructed 
array. 



Proof: These properties follows immediately from the 
fact that the entries of S by the order of the folded-row are 
consecutive elements of the M-sequence S. The two cyclic 
shifts of S have the same folded-row up to a cyclic shift. 
Therefore, these four properties are a direct consequence 
from the related properties of the M-sequence. ■ 
Lemma 21: Let A be a lattice tiling for the shape S and 8 
be a direction for which the triple (A, S, 8) defines a folding. 
Let B be a binary sequence of length |«S|. Let Pi and P2 be 
two points for which Pi — c(Pi) — P2 — c(P2). Then, for 
any two positive integers fci and k 2 the two fci x &2 windows 
of (A, S, 8, B) whose leftmost bottom points are Pi and P2 
are equal. 

Proof: The lemma is an immediate consequence from 
the definition of the lattice coloring induced by (A, S, 8) and 
the definition of (A, S, 8,B). ■ 

Theorem 16: Assume A define a lattice tiling for an n\ x 
712 array A, such that 711712 = 2 klk2 — 1. Assume further that 
A defines a lattice tiling for the shape S and (A, S, 8) defines 
a folding for the direction <5. Then, if we fold an M-sequence 
S into S in the direction 8, the resulting shape S has the 
ki x fc 2 window property if and only if the 7i! x n 2 array 
A has the fci x fc 2 window property by folding S into A in 
the direction 8. 

Proof: Since A is a lattice tiling for both A and S 
there is a sequence of arrays A® = A, Ai,...,A r — S, such 
that \A l+1 \Ai\ = \Ai \ A i+1 \ = 1, < i < r - 1, A is a 
lattice tiling for Ai, < i < r, and the origin is contained in 
Ai, < i < r. Moreover, it is easy to verify that given the 
shape Ai, Pi = Ai+i \ Ai, P2 = Ai \ Ai+\, we have that 
P2 = P\ — c(P{) with respect to Ai. The theorem follows 
now by induction and using Lemma |2"T1 ■ 

Theorem [16] does not give any new information about 
window sizes which are not covered in [40], [41]. The 
following lemma provides such information. We say that 
a shape S of size 2™ — 1 has the Q window property if 
|<2| = n and each nonzero value for Q appears exactly once 
in a copy of S, where S is considered to be a cyclic shape. 

Lemma 22: Let A be a lattice tiling for a shape S, \S\ = 
2™ — 1, 6 be a direction vector, and S be an M-sequence of 
length 2" — 1. Let Q be a shape with volume n. If in the 
array S' defined by (A, S, 8, S) there is no copy of Q which 
contains only zeroes then S has the Q window property. 

Proof: By the Shift-and-Add property, S' has two 
identical copies of Q if and only if S' has a copy of Q 
which contains only zeroes. Thus, S' has the Q window 
property if and only if there is no copy of Q in S' which 
contains only zeroes. ■ 
We can use now the properties we have found for the 
generalized folding to obtain various results. An example 
is given in the following corollary. 

Corollary 10: Let A be a lattice tiling for a shape S, 
\S\ = 2™ - 1, and S be an M-sequence of length 2" - 1. If 
2™ — 1 is a Mersenne prime then (A, S, 8, S) has the 1 x n 
and the n x 1 window property for any given direction vector 
8. ' 

Example 5: Consider the following M-sequence S = 
0000100101100111110001101110101 of length 31. Let A 
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be a lattice tiling for a corner CR(5,7;1,4) with the generator 
matrix 

3 4 
10 3 



G 2 



By folding of S in the direction (+1,0) we obtain the 
following pseudo-random array 



1 





1 




1 


1 





1 


1 


1 





1 


1 


1 


1 











1 





1 


1 








1 














1 









This array has the 5 x 1 and 1x5 window properties. 
Out of the 19 shapes of size 5 with exactly two rows it does 
not have the window property only for the following three 
shapes: 



d c 



i 



The pseudo-random array obtained by folding S by the 
direction (0, +1) is 



1 


1 


1 










1 


1 





1 


1 














1 


1 


1 





1 





1 





1 


1 





1 











1 






It has the 5x1 and 1x5 window properties. But, out of 
the 19 shapes of size 5 with exactly two rows it does not 
have the window property for eight shapes. 

Both pseudo-random arrays have a window property for 
the star shape given by 



IX. Conclusion and Open Problems 

The well-known definition of folding was generalized. 
The generalization and its applications led to several new 
results summarized as follows: 

1) The generalization is based on a lattice tiling for a 
shape S and a direction 5. The number of possible 
nonequivalent directions is , Necessary and suf- 
ficient conditions that a direction defines a folding are 
derived. 

2) Folding a i?2-sequence into a shape S result in a 
distinct difference configuration with the shape S. 

3) Lower bounds on the number of dots in a distinct 
difference configuration with shape of regular poly- 
gon, circle, and other interesting geometrical shapes 
are derived. 

4) Low redundancy multidimensional codes for correct- 
ing a burst of length two are obtained. 



5) New pseudo-random arrays with window and corre- 
lation properties are derived. These arrays differ from 
known arrays either in their shape or the shape of their 
window property. 
The discussion on these results leads to many new inter- 
esting open problems. We conclude with a list of six open 
problems related to our discussion. 

1) We have discussed several applications for the folding 
operation in general and for the new generalization 
of folding in particular. We believe that there are 
more interesting applications for this operation and we 
would like to see them explored. 

2) The construction for DDCs whose shape is a quasi- 
perfect hexagon works for infinite number of param- 
eters. But, the set of parameters is very sparse. Its 
density depends on the number of primes obtained 
by Dirichlet's Theorem. This immediately implies the 
same for the parameters of DDCs whose shape is a 
regular polygon. We would like to see a construction 
of such DDCs with a dense set of parameters. 

3) What is the lower bound on the number of dots in 
a DDC whose shape is a circle with radius Rl We 
conjecture that the lower bound is -JtxR + o(R). 

4) We would like to see an asymptotic improvement on 
the lower bounds on the number of dots in a DDC 
whose shape is a regular n-gon with radius R. 

5) Are there cases where we can improve the upper 
bound on the number of dots in these DDCs asymp- 
totically? 

6) We would like to see a more general theorem which 
connects folding of M-sequences and general window 
property. 

Appendix A 

In this Appendix we prove the necessary and sufficient 
condition for a triple (A, S, 5) to define a folding. For the 
proof of the theorem we use the well known Cramer's 
rule [61] which is given first. 

Theorem 17: Given the following system with the n 
linear equations and the variables xi, X2, ■ ■ ■ , x n 



an 


au ■ 


■ ain 




Xi 




' h ' 


an 


a 2 2 ■ 


■ a2n 




x 2 




b 2 


a n i 


a n 2 


ann 











If 



A = dct 



an ai2 
«2i a 22 

a n i a n2 



air. 

a2r. 



then Xk — fov 1 < k < n, where 



det 



an • ■ • ai(fc-i) bi Oi( fe+ i) . . . 

«21 • ■ • 02(fe-l) b 2 «2(fc+l) • ■ ■ 

a„i . . . a„(fe_i) b n a n (k+i) 



air. 

a2r. 
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Let A be a _D-dimensional lattice tiling for the shape S. 
Let G be the following generator matrix of A: 



G = 



"11 "12 
«21 "22 

"£>1 "D2 



"l£> 
"2D 

VDD 



Given the direction vector J = (di, tfe, • • • , do), w.l.o.g. 
we assume that the first l\ >\ values of 6 are positives, 
the next £2 values are negatives, and the last D — £\ — £2 
values are O's. By Lemma [5] and Corollary if (A, S, S) 
defines a folding then there exist D integer coefficients 
ai, as, ... , a D such that 

D 

^aj{v j i,v j 2,...,Vj D ) = 
3=1 

(\S\di, . . . , |<S|d^,-|<S|d^ + i, . . .,-\S\d il+i2 ,0, ... ,0), 

and there is no integer i, < i < \S\, and D integer 
coefficients (3i , (32 , ■ ■ ■ , (3d sucn that 



3=1 



(i ■ di, . . . ,%■ di x ,— i- di x . 



-i ■ di 



,0) 



Hence we have the following D equations: 

D 



3=1 



\S\-d r , l<r<4, 



(2) 



D 
3=1 



\S\-dr, ei + l<r<h+h, (3) 



Y ctjVjr =0, £1 + £2 + 1 < r < D . (4) 
3=1 

Let r = di if £1 +^2 = 1 and r = g.c.d.((ii, dsi ■ ■ ■ , <^i+£ 2 ) 
if £\ +£2 > 1. The I? equations in ©, (0, (|4]l are equivalent 
to the following D equations: 

D 

^ajVji = \S\ - di, 

3=1 



3=1 



0, 2<r<4, 



D 

£ 

3=1 



Q' , 



d\Vj r + d r Vji 



0, £1 + 1 < r < £1 + £ 2 



a^ ir = 0, £1 + £ 2 + 1 < r < D. 

3=1 

We define now a set of D(D — 1) new coefficients u r j, 
2 < r < D, 1 < j < D, as follows: 



Urj = 1 Jr ^ J1 for £1 + 1 < r < ti + £2, 



u r j = Vj r for l\ + £ 2 + 1 < r < D. 
Consider the (D — 1) X D matrix 



H 



U21 "22 

"31 "32 
UD1 U D 2 



U 2 D 
"3D 

UDD 



Using Theorem[T7]it is easy to verify that the unique solution 
for the afc's is 



Oik 



-1) 



^diT^+^-Metiffc 



4 



l+ia-l 



(5) 



where Hk is the (£) — 1) X (D — 1) matrix obtained from 
H by deleting column k of H. 

Lemma 23: For each k, 1 < fc < D, t divides ctk defined 
in ©. 

Proof: Consider the following D x D matrix 



G 



Vll 


"21 ■ 


■ VD1 


U21 


W22 ■ 


■ U 2 D 


M31 


"32 ■ 


■ u 3D 



UD1 U D 2 



UDD 



By the definition of the entries in the matrix H and since 
detG = \S\ it follows that that dctG = \S\ . 
detG in Theorem [TTI is equal A, while Ak is equal \S\ • 
di (—) 1 2 Y, for some integer Y. Therefore, = tY 
and the lemma follows. ■ 

This analysis leads to the following theorem. 

Theorem 18: If A is a lattice tiling for the shape S 
then the triple (A, S, S) defines a folding if and only if 
g.c.d.(^, a .. . , aa) = 1 and g.c.d.(r, \S\) = 1. 

Proof: Assume first that (A, S, 8) defines a folding. 

Now, assume for the contrary that 
g.c.d.(^, J* .. . , aa) = n > 1 or g.c.d.(r, = ^2 > 1. 
We distinguish between two cases. 
Case 1: 

Assume that g.c.d.(^, sa, . . . , ^a) = i/ x > 1. 

Equations (fJJ, OJ, and © have exactly one solution for 
the ai's given in ©. Since g.c.d.(^, Sa, . . . , ^a) = Vlj j t 
follows that (3i — 1 < i < D, are integers. Therefore, 
we have 



151 

yBjVjr = —d r , l<r<£i, 

3 = 1 



D 

£/^. 
3=1 



,, - — —dr, £1 + 1 < r < 4 + £ 2 , 

TVl 



D 



Pjv jr = 0, ii + £ 2 + l<r<D, 

3=1 
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i.e. 



, ° \S\, \S\ \S\ 

TV\ TV\ TV\ TV\ 

and as a consequence by Lemma [3] we have that (A, S, S) 
does not define a folding, a contradiction. 
Case 2: 

Assume that g.c.d.(T, |<S|) = zaj > 1. 
Let $ = gi, 1 < i < £i + 1 2 . Therefore, 

D 

( — di, . . . , — dt x , <k 1 +i, ■ ■ ■ , ae 1+ e 2 , u, . . . , U), 

Z/ 2 1^2 ^2 ^2 

and as a consequence by Lemma [5] we have that (A,S,S) 
does not define a folding, a contradiction. 

As a consequence of Case 1 and Case 2 we have that 
if (A, S, 8) defines a folding with the ternary vector 6 then 
g.c.d.(^, sa . . ., aa) = 1 and g.c.d.(r, |5|) - 1. 

Now assume that g.c.d.(2i, 2a, . . . , 22.) = 1 and 
g.c.d.(r, |«S|) = 1. Consider the set of Z? equations defined 
by 

D 

^ Qij (Uj-i ,V j2 ,..., Vjp) = (6) 
3=1 

(\S\di, \S\d ei ,-\S\de 1+ i, . . . , -|<S|^ 1+ ^,0, ... ,0), 

Since the rows of G are linearly independent, it follows that 
this set of equations has a unique solution for the ctj's (but, 
these coefficients are not necessary integers). Using the same 
analysis proceeding the theorem, we have by the Cramer's 
rule that this solution is given by (0) and hence the aj's 
are integers. Assume for the contrary that (A, S, S) does not 
define a folding. Then, by Lemma|5]we have that there exist 
D integers /3,-, 1 < i < D, such that 



D 

^2Pj(Vjl,Vj2,..-,V jD )= (7) 

i=i 

(£■ di, ...,£■ d£ 1} -£■ di 1+ i, ...,-£■ di 1+ e 2 ,0, ... ,0), 

for some integer < £ < \S\. 

Since the rows of G are linearly independent then there 
exists exactly one set of /3j's (integers or non-integers) which 
satisfies Q. Let v — g.c.d.(£, \S\), where clearly 1 < v < 
£ < \S\. From equations © and (|7]) we obtain 

D 

^r(£a j )(v j i,v j 2,...,v jD ) = 

3=1 

(^|5|di, . . . , l\S\d tl ,-l\S\dt 1+u -l\S\d tl+ta ,0, . . . , 0) 

D 

= ^2(\S\0j)(Vjl,Vj2, . ■ ■ ,V jD ) , 
3=1 



' Fig. 7. A comer CR(7, 11; 2, 4) 



Since the rows of G are linearly independent it implies 
that £a t = \S\f3, for each 1 < i < D, i.e., ft = fe. 



|5| 

ijv 

\S\I 



f is an integer and v — g.c.d.(£, |<S|) implies that 
-O-i, 1 < i < D. g.c.d.(£/v, \S\/v) = 1 and hence 



ft 

— divides a$ for each £, 1 < i < D. g.c.d.(r, |5|) = 
1, r divides a,, and hence — divides — for each L 
1 < i < D. Hence, g.c.d.(^, 2a, . . . , 22) > Jfl. But, 
g.c.d.(2i,^,...,^) = 1 and hence =]<S|, i.e., £ > 
a contradiction. Thus, (A, S, S) defines a folding. ■ 

Appendix B 

In this appendix we consider DDCs with two special 
shapes, called corner and flipped T. The DDCs with these 
shapes and special parameters are important in applying 
Theorem [6] to obtain other DDCs such as triangles in the 
square grid and hexagonal spheres in the hexagonal grid. 

A. Corner 

A corner, CR(hi -\-h%, w± + w 2 ; h,2, W2), is an (hi + /12) x 
(wi + W2) rectangle from which an /12 X W2 rectangle was 
removed from its right upper corner. An example is given 
in Figure [7] Let S be a CR(h\ + h 2 , wi + W2] h 2 , W2) and 
let A the lattice with the following generator matrix 



G = 



wi hi 
—W2 hi + h 2 



Clearly, A is a lattice tiling for S. A general result 
concerning DDCs whose shape is a corner seems to be 
quite difficult. We will consider the case which seems to 
be the most useful for our purpose. First note, that by 
Theorem [TJ S = (0, +1) defines a folding for A if and 
only if g.c.d.(uii, w 2 ) = L Assume first that hi = h 2 
and \wi — w 2 \ < 3. By Theorem [TT] we have an ni x n 2 
rectangle Q such that nin2 = p 2 — 1 for some prime p, 
2ni ^ fei+ft2 an( j n j s even> Now, we will make new 
choices for hi, h2, w%, and W2, which are close to the old 
ones. Let hi = /12 = nil we distinguish between three cases 
of n 2 : 

(W.l) If n 2 = 3w + 1 then 101 = oj and W2 = w + L 
(W.2) If n 2 = 3u> + 2 then Wi = uj + 1 and W2 = w. 
(W.3) If n 2 = 3w then we distinguish between two cases: 

• if uj — 1 = (mod 3) then w\ = uj + 1 and W2 = 
uj-2. 

• if uj — 1 ^ (mod 3) then wi = uj — 1 and W2 = 
uj + 2. 

It is easy to verify that the size of the new corner CR(/ii + 
h 2 ,wi + w 2 ;h 2 , W2), S', is ni7i2 = p 2 — 1, A is a lattice 
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Fig. 8. A flipped T FT(5, 17; 4, 6) 



tiling for S', (A, S', 5), 5 = (0, +1), defines a folding, and 
we can form a doubly periodic 6>'-DDC with it. Hence, we 
have the following theorem. 

Theorem 19: Let n\ and be two integers such that 
n\rii = p 2 — 1 for some prime number p, ri2 = 2u>i + W2, 
where ri\ is an even integer, w\, W2, are defined by (W.l), 
(W.2), (W.3). Then there exists a doubly periodic 5-DDC, 
whose shape is a corner, CR(2ni, w\ + W2; 711,1112), with p 
dots. 



B. Flipped T 

A flipped T, ¥Y{h, wi + w 2 + w 3 ;wi,W3), is an (2h) x 
(■W1+W2+W3) rectangle from which an/ixwi rectangle was 
removed from its left upper corner and an h x W3 rectangle 
was removed from its right upper corner. An example is 
given in Figure [8] Let S be a FT(h, w\+W2+ W3; w\, W3) 
and let A the lattice with the following generator matrix 



G = 



W\+ W2 h 

wi + 2u>2 + W3 



Clearly, A is a lattice tiling for S. A general result 
concerning DDCs whose shape is a flipped T seems to 
be quite difficult. We will consider the case which seems 
to be the most useful for our purpose. First note, that by 
Theorem Q] 5 = (0, +1) defines a folding for A if and only 
if g.c.d.(u>i +W2, wi + 2w2 + W3) = 1 which is equivalent to 
g.c.d.(u>i + W2, W2 + W3) = 1. Assume that \w± — ws\ < 4. 
By Theorem [TT] we have an n\ x 122 rectangle Q such that 
nin2 = p 2 — 1 for some prime p, — ss — . „ h , — , and no 
is even. Now, we will make new choices for h, w\, and W3, 
which are close to the old ones. Let /1 = ni; we distinguish 
between two cases of n.2'- 

(Y.l) If ri2 = 4w then w\ — 2u + 1 — W2 and 103 = 2u ~ 

1 - U>2- 

(Y.2) If ri2 = iuj + 2 then w\ = 2lu + 3 — W2 and W3 = 

2iU — 1 — W2- 

It is easy to verify that the size of the new flipped T, 
FT(/i, w\ + W2 + W3; w\, W3), S', is ni?i2 = p 2 — 1, A is a 
lattice tiling for S' , (A, 5', 8), 5 = (0, +1), defines a folding, 
and we can form a doubly periodic <S'-DDC with it. Hence, 
we have the following theorem. 

Theorem 20: Let n\ and W2 be two integers such that 
ri2 = W\ + 2u>2 + W3, w\, W3, are defined by (Y.l), (Y.2), 
and n\n2 = p 2 — 1 for some prime number p. Then there 
exists a doubly periodic 5-DDC, whose shape is a flipped 
T, FT(ni, wi + u>2 + W3; W\, W3), with p dots. 



Appendix C 

In this section we demonstrate how Theorem [6] is applied 
for several geometric shapes (having the role of Q in the 
theorem), where our shape S in the doubly periodic 5-DDC 
is an appropriate corner, flipped T, or quasi-regular hexagon. 

C. Equilateral Triangle 

Let Q be an equilateral triangle with sides of length B. 
The area of Q is ^B 2 and hence an upper bound on the 

number of dots in Q is ^-B + o(B) = 0.658B + o(B). For 

our shape S we take a flipped T, FT^, ^§5; ^) 
which overlaps in its shorter base with the base of Q. These 
bases of S and Q share the same center. The area of S 
is ^B 2 and hence the density of the array is — ^— . The 

4 34 B 

intersection of S and Q, A(Q,5), equal to hH^llB 2 . 
Therefore, a lower bound on the number of dots in Q is 
3V2-2V3 B + ( B \ = 0.59165 + o(B) and the resulting 

packing ratio is 0.899 . The same result can be obtained by 
using other structures instead of a flipped T. 

D. Isosceles Right Triangle 

Let Q be an equilateral triangle with base and height of 
length B. The area of Q is ^B 2 and hence an upper bound 
on the number of dots in Q is ^B+o(B) = 0.707B+o{B). 

For our shape S we take a corner CR(^J~^B, \J~^B; -^=, ^=) 
which overlaps in its two shorter sides with the base and 
height of Q. S and Q shares the intersection vertex of these 
sides. The area of S is ^B 2 and hence the density of the 
array is The intersection of S and Q, A(Q,S), equal 
to (y/6 — 2)B 2 . Therefore, a lower bound on the number of 
dots in Q is (Vl2 - 2V2)B + o(B) = 0.63567B + o(B) 
and the resulting packing ratio is 0.899 (exactly as in the 
case of an equilateral triangle). 

E. Regular Pentagon 

Let Q be a pentagon with radius R. The area of Q is 
I sin 4? and hence an upper bound on the number of dots 
in Q is 1.54196i? + o(R). Let S be a quasi-perfect hexagon 
having a joint base with Q and two short overlapping 
sides with Q, where these sides are connected to this base 
(see Figure |5}. The distance between the base and the 
diameter of S is aR, 2 sin cos ^ < a < (1 + sin )/2. 
The length of the base is 2R sin ^ and the length of the 
diameter of S is 2i?sin^ + 2ai?tany^. Hence, the area 
of S is ^sin^ + 2a tan j^)aR 2 and the density of the 



array is 



The area of the intersection 



^/ 4a sin ^ +2a 2 tan H 

between Q and S, A(S, Q), is computed by subtracting 
from the area of S the area of the two isosceles triangles 
<7i and (T2. The lower bound on the number of dots is 
, 1 „ A(6>. Q). The maximum on this lower 

^/4asin f +2a 2 tan j^R y ' > 

bound is obtained for a = 0.814853, i.e., the lower bound 
on the number of dots in a pentagon with radius R is 
1.45992i? + o{R) yielding a packing ratio of 0.946795. 
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aR < 



Fig. 9. Quasi-regular hexagon intersecting a regular pentagon 

F. Regular Heptagon 

Let Q be a regular heptagon with radius R. The area of 
Q is | sin ^-R 2 and hence an upper bound on the number 
of dots in Q is 1.65421i? + o(R). Let S be a quasi-perfect 
hexagon constructed as follows. We refer to the sides of 
Q as side 0, side 1, side 2, side 3, side 4, side 5, side 6, 
in consecutive order clockwise. Let's denote the six sides 
of S by side A, side B, side C, side D, side E, side F, 
in consecutive order clockwise, where side A is the lower 
base of S. Sides B and C of S overlap sides 1 and 2 of Q, 
respectively; sides B and C are longer than sides 1 and 2. 
The two bases of S, sides A and D, have angles with 
sides B and C, respectively. Side A intersect sides and 6 
of Q; side D intersect sides 3 and 4 of Q. The length of 
the segment, on side 0, from the vertex of the intersection 
between sides and 1 and the intersection of side A and side 
is xR. Finally, sides E and F of S are parallel to sides B 
and C, respectively; Side E intersect sides 4 and 5; side F 
intersect sides 5 and 6. The distance between the vertex of 
the intersection between side E and F of S and side 5 of 
Q is aR. Computing |<S|, A(<S, Q), and the lower bound on 
the number of dots in Q, i.e., ^Sl, as functions of x and 



a implies that the maximum is obtained for x — 0.432042 
and a = 0.0840633, and the lower bound on the number of 
dots in Q is 1.58844i? + o(R) yielding a packing ratio of 
0.960241. 

G. Regular Octagon 

Let Q be a regular octagon with radius R. The area of Q is 
4 sin jR 2 and hence an upper bound on the number of dots 
in Q is 1. 68179-R + o(R). Let S be a quasi-perfect hexagon 
having a joint diameter of length 2R with Q and overlapping 
four side edges with the Q. The distance between the diame- 
ter of the hexagon (octagon) and a base of S is aR. The area 
of the S is 4aR 2 — 2a 2 ^"j^ i? 2 and hence the density of the 



array is 



4a-2a 2 "".7, R 



The intersection between Q and 



S, A(S, Q), is 4 sin \R 2 - 2(1 - a) 2 R 2 



Therefore, 



a lower bound on the number of dots in the octagon is 

4sin f fl-2(l-a) 2 _R 5 & '°T 

. . The maximum is obtained for a = 



ia~2a 2 



1 ~5~ 



0.872852 and hence a lower bound on the number of dots 
is 1.62625i? + o(R) and the packing ratio is 0.966977 . 



H. Regular Nonagon 

Let Q be a regular nonagon with radius R. Let S be a 

BL ft 

77 -it, 



quasi-regular hexagon with radius p, where p = - 
such that Q and S share the same center and there is an 
overlap in three pairs of edges between Q and S. The 
area of Q is | sin ^-R 2 and hence an upper bound on the 
number of dots in Q is 1.700748i? + o(R). The area of 

S is ( S 'sTn? ) ^ 2 anc * hence the density of the array 

is — [ — . The area of the intersection between Q 

3^\/3,Rsin ^§- 

and S, A(S, Q), is ( s -^) ' R 2 - 6 s -^^- R 2 = 

A \ Sm ~3 I Sln ~3 

2.8625667i? 2 . Therefore, a lower bound on the number of 
dots in the nonagon is 1.63672i?+ o(R) and the packing 
ratio is 0.96235 . 



/. Regular Decagon 

Let Q be a regular decagon with radius R with sides 0, 1, 
2, 3, 4, 5, 6, 7, 8, 9 in consecutive order clockwise. The area 
of Q is 5 sin ^R 2 and hence an upper bound on the number 
of dots in Q is 1.71433i? + o(R). Let S a quasi-perfect 
hexagon with sides A, B, C, D, E, and F, in consecutive 
order clockwise, where A is the lower base of S. Sides B 
and C of S overlap with sides 1 and 3 of Q; sides E and 
F of Q overlap with sides 6 and 8 of Q. The two bases A 
and D of S have distance aR to the diameter of S which 
connects the intersection vertex of sides B and C with the 
intersection vertex of sides E and F. The distance between 
the diameter and a base (A or D) is aR. The area of S is 
s = 2(2 sin 4f + 2 si "^^" f - ^4ra)aR 2 and the density 

sm — sin — ^ ^ 

of the array is 4-. Finally, A(<S, Q) = (5 sin \ -2S^(1- 
a) 2 )R 2 . The lower bound of the number of dots in Q is 
-4= A(<S, Q). The maximum on this lower bound is obtained 
for a = 0.923286; the lower bound is 1.64786.R + o(R) and 
the packing ratio is 0.961229 . 
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